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INTRODUCTION 

Depuis plusieurs années, on s’efforce d'utiliser le plus possible des 
méthodes mathématiques dans la manipulation des données biologiques. 
Ces méthodes sont principalement employées dans l’étude des populations 
animales ou végétales et dans leurs relations avec l'environnement. 

Une première approche consiste à essayer d’inférer une équation mathé¬ 
matique à partir d’une courbe expérimentale, tracée à partir de résultats 
biologiques. Une équation est ainsi élaborée pour que sa représentation 
graphique soit comprise dans la surface globalisant les données de l’expé¬ 
rience. Un test de corrélation effectué permet rarement de constater une 
bonne correspondance entre les deux tracés. Bien mieux, habituellement, 
il est possible de construire un nombre quelconque de courbes différentes 
rendant compte de façon approximative des valeurs des résultats expéri¬ 
mentaux. Il est clair qu’une telle équation peut être considérée, au maxi¬ 
mum, comme une moyenne des données des expériences. Elle ne peut, 
en aucun cas, représenter une « loi » biologique et les différents para¬ 
mètres mathématiques utilisés n’ont en soit aucune signification biologique 
(voir, par exemple, les travaux de : Burnett, 1953b; Leon, 1968; Utida, 
1941a, b. c, d, e, 1942a, b, 1943a, b, c, 1950). 

Heureusement, ces applications naïves des mathématiques à la biologie 
ne sont pas les seules tentatives effectuées pour quantifier davantage les 
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observations et les concepts biologiques. En effet, une autre voie existe. 
Elle consiste en l'élaboration de modèles. 

Dans ce cas, des facteurs biologiques simples, les composants de base 
d'un phénomène, sont représentés par une série de symboles. Ces facteurs 
sont associés dans une équation où dans une série d'équations exprimant 
des relations biologiques. 

Par des manipulations et des additions appropriées, on obtient ainsi une 
représentation mathématique, un modèle, équivalente aux phénomènes 
biologiques à étudier. Dans ce cas, et ceci est très important, les para¬ 
mètres des équations ne sont pas quelconque mais représentent des fac¬ 
teurs biologiques ou leurs relations. 

L'utilisation de ces équations permet de prédire des effets suivant 
certaines causes, même en dehors des limites des données biologiques 
généralement disponibles. Les extrapolations effectuées permettent, d'une 
part, de vérifier les hypothèses de base et, d'autre part, de mettre en évi¬ 
dence de nouveaux phénomènes. 

Un certain nombre d'équations de ce genre ont été construites. Nous 
allons étudier systématiquement les principaux modèles proposés et voir 
sur quelles bases biologiques ils reposent. Nous partirons des plus simples 
pour en arriver aux plus complexes. 

Nous envisagerons ainsi successivement : 

I. Les modèles déterministes. 

1. Population croissant indéfiniment. 

2. Population affectée d’un coefficient de natalité et de mortalité. 

3. Population croissant dans un espace limité. Equation de Verhulst. 

4. Croissance de deux populations dans un espace limité. 

5. Relation entre proies et prédateurs. Equations de Lotka-Volterra. 

6. Modèles à multi-relations. Méthode de Garfinkel ef al. 

II. Les modèles stochastiques. 

1. Population croissant dans un espace limité. Equation modifiée de 
Verhulst. 

2. Relations entre proies et prédateurs. Equations modifiées de Lotka- 
Volterra. 

Le modèle que nous proposons plus loin se différentie des tentatives 
précédentes par un parti pris, à la base, de simplicité. Les équations que 
nous proposons ne sont pas immédiatement déduites de résultats expéri¬ 
mentaux. Elles ne sont pas non plus construites à partir de globalisations 
de multiples paramètres (coefficient de croissance, par exemple). Au 
contraire, notre modèle est constitué d’un grand nombre d'équations 
réduites représentant chacune un phénomène biologique évident. Par 
exemple; il est bien certain que les animaux se déplacent à une certaine 
vitesse V. 
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Cette succession d'équations simples est alors intégrée dans un modèle 
global où 1 on calcule, à l'aide d’un ordinateur, les relations entre les 
divers paramètres. Les résultats peuvent sortir suivant diverses formes dont 
la plus commode semble être la représentation graphique des nombres 
d'animaux de chaque population, en fonction du temps. 

Le but de l'étude est de posséder ainsi un modèle général ou tous les 
paramètres sont parfaitement caractérisés. Puis d'expérimenter sur ce 
modèle pour mettre en évidence les interactions lointaines entre les 
variables. 

Le travail présenté ici détaille les paramètres et les méthodes utilisées. 


I. LES MODELES DETERMINISTES 

1. Population croissant indéfiniment 

Les propositions utilisées dans ce modèle très simplifié sont les sui¬ 
vantes (Yule, 1924) : 

— les organismes sont considérés comme étant immortels; 

— leur coefficient de reproduction est le même pour chaque individu 
et ne varie pas avec le temps; 

— il n'y a pas d’interaction entre les individus; 

— l'espace dans lequel les animaux peuvent vivre n’a pas de limite. 

Si N, = nombre d’organismes au temps t 
À = coefficient de reproduction 

On a : 

c/N, 

-- AN, 

dt 

Si la population au temps O était N,„ on en déduit : 

N, = N 0 eX‘ 


La croissance de la population est donc exponentielle, dans ces condi¬ 
tions. 

Au point de vue biologique, cette équation a évidemment une portée 
très limitée car les propositions de base sont trop simples. Elle n'est 
d'application — limitée dans le temps — que dans le cas de certaines 
croissances cellulaires. 

2. Population affectée d'un coefficient de natalité 

ET DE MORTALITÉ 

Les propositions utilisées sont les mêmes que dans le modèle précédent 
avec l'addition d'un coefficient supplémentaire de mortalité n. 
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On a immédiatement : 

N, = N„ e(A — n)t 

L’utilisation de cette équation est aussi limitée que la précédente. 


3. Population croissant dans un espace limité. 

Equation de Verhulst 

Ce modèle est le plus ancien qui ait jamais été présenté puisqu'il est 
dérivé des travaux de Verhulst (1838). 

Dans ce modèle, les organismes sont considérés comme ne pouvant 
croître indéfiniment : leur coefficient de reproduction est fonction de 
la densité de la population : A = / (N). 


On a : 


d N 

-= N / (N) 

dt 

Si l’on considère que A est une fonction linéaire de N, on peut écrire : 


D'où : 


A = f (N) = a — h N avec a et b > O 


dN 

-= N (a - bN) 

dt 

Qui peut s’écrire : 

dN b dN 

-h- = a dt 

N (a - bN) 


(3.1) 


En intégrant, on obtient : 

N, 

-= C e‘" 

(a - b N,) 


C étant une constante d’intégration. 

D'où : 

C a e'" a/b 

N, =--- 

1 + C b e'" 1 + e-'-'/Cb 

En choisissant C tel que Cb = e (Bartlett, 1960), on a : 

a/b 


N, 


1 + e-“ l '-V 
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La valeur N, tend donc vers une limite — a/b ■—, quand t tend vers 
l'infini : la population atteint un niveau de saturation. 

Cette équation a été utilisée de nombreuses fois car, dans certains cas, 
il existe une bonne concordance entre les résultats expérimentaux obtenus 
et le modèle (Gause, 1934; Lotka. 1925; Odum, 1959; Pearl, 1924, 
1932, 1934; Pearl et al.. 1920, 1922). 

Examinons, maintenant, le postulat biologique de base utilisé dans ce 
modèle : le coefficient de reproduction d'une population dépend, d’une 
façon linéaire, de la densité des organismes. 

Ce postulat pèche malheureusement par une trop grande simplicité. 
Les recherches de Stanley (1932 a et b, 1934) montre ainsi que le 
nombre d’œufs pondus par Tribolitim confusum Duv., dans des conditions 
expérimentales contrôlées, dépend d'au moins 14 facteurs. Les ramener 
à un seul est évidemment une utopie. 

Néanmoins, de par sa flexibilité — il suffit de donner des valeurs 
appropriées à a et b ou d'ajouter des constantes, l'équation de Verhulst 
a été largement utilisée. Ullyett (1953) a discuté longuement les nom¬ 
breux manques de ce modèle qui peut être considéré, en fait, comme une 
simple formulation mathématique d’une courbe empirique, les facteurs 
utilisés dans les équations n'ayant une signification biologique imprécise, 
malgré des formulations nouvelles telles celles de Lewis (1942), Lefko- 
vitch (1965, 1966). Leslie (1945. 1948, 1957, 1958, 1959, 1960), 
Leslie et al. (1940, 1949). Lewis (1942), Pollard (1966) et Vander- 
meer (1969). 


4. Croissance de deux populations dans un espace limité 

Ce modèle est directement déduit de l'équation précédente (§3) mais 
sa formulation en est plus complexe. On va considérer, en effet, que le 
coefficient de reproduction d’une espèce 1 est limité par la densité de la 
population 1 mais aussi par la population 2. 

On a les paramètres suivants : 

espèce 1 espèce 2 

b iX 
b„ 

b,r 

a, b u , a-, et b 22 sont semblables à a et b du § 3. 

by. exprime la façon dont la présence de l'espèce 2 affecte la croissance 
de l'espèce 1. 

b 21 exprime la façon dont la présence de l'espèce 1 affecte la croissance 
de l’espèce 2. 
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On obtient ainsi les 2 équations différentielles suivantes, Ni et N 2 
étant le nombre d'individus de l'espèce 1 et de l'espèce 2. 

cfN, 

-= N, (a, - b u N, - b,. NJ 

dt 

rfN. 

-= N, (a 3 - b 2i N, - b,, NJ 

dt 

Ces deux équations n'ont pas de solutions immédiates. On arrive 
cependant à les résoudre en introduisant une simplification supplémen¬ 
taire. On va considérer que l'effet inhibiteur sur la croissance de la 
densité totale des populations 1 et 2 est identique sur 1 et 2. 
En d’autres termes, chaque individu de l'espèce 1 ou 2 entre en compé¬ 
tition avec une population de dimension N = N, + pN 2 , p étant un fac¬ 
teur indiquant que les deux espèces peuvent avoir un effet inhibiteur 
différent. 

Comme dans le cas 3.1. le système d'équations différentielles peut 
alors s'écrire : 


d N, 

- N, (a, 
dt 

- b t N) 

(4.1) 

</N.. 

-= N., (a.. 

dt 

- bi N) 

(4.21 


Ce qui revient à dire que : 

bu N, + b v . N. .. = b, N 
b,, N, + b„ N, = b, N 

En remplaçant N par sa valeur, on a : 

b u N, + bu N s = b, (N, + pNJ 
bu N, + b y, N» = bi (N, + pNJ 

D'où : 

N, (bu - b ,) + N,, (b,, - pb t ) = 0 
N, {bu — ij + N,, (bu - b..p) = 0 

Comme N, et N 2 sont positifs, on a : 

£»i i = b , 
bu = bi 

bu bu 

-= p ou -= p 

b, bu 
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b.,; b,-, 

- = p ou -= p 


b 2l — b ii f > 22 

Les équations 4.1 et 4.2 peuvent alors s'écrire : 
c/N, 1 d 

-.-= — ln N, = a, — b, N 

dt N, dt 

c/N 2 1 d 

-.-= — lnN. = a 2 — b-. N 

dt N 2 dt 


En éliminant N. on a : 

d c/N, 

— (b-, l/i N, — b, ln N 2 ) = — ln (-) = a, b-, — a 2 />, 

dt dt ’>i 

N. 2 

Si N,, 1 et N,, 2 sont les populations de l'espèce 1 et de l'espèce 2 au 
temps O. on a : 

N, (N,, 1 ) (a, 6, — a 2 b,)/ 



Si a, b-, > a 2 b , N 2 tend vers O 

Si a, b-. < a-, b, N, tend vers O 

La restriction originelle b v .b n = b n b,, peut être modifiée et devenir 

b la b-n b i, b 2i . 

Dans ces conditions, 4 cas peuvent se produire suivant les valeurs 
choisies (Slobodkin, 1961) : 

— l’espèce 1 disparaît; 

— l’espèce 2 disparaît; 

— il apparaît un équilibre stable; 

— il apparaît un équilibre instable. 

Ces équations et leurs résultats peuvent apparaître en biologie car 
elles semblent, par la variété de leurs résultats, représenter davantage 
les variations des populations d'organismes. Mais il s’agit d'une illusion. 
Les coefficients a, b, b„, etc, n'ont qu'une signification biologique très 
vague et regroupent un grand nombre de phénomènes qui ne sont 
pas constants. 
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Théoriquement, cependant, ces équations pourraient être de quelque 
utilité dans le cas de populations expérimentales croissant dans des 
conditions d'une grande stabilité. C'est ce que Gause (1935) et Gause 
et al. (1935) avaient obtenu avec des Protozoaires. Cependant, dans le 
cas d’insectes, il faut vraiment être de bonne volonté pour voir une 
relation bien définie entre les résultats et les prévisions des équations 
(Crombie, 1945, 1946). 

Dans le cas des Poissons. Walters (1969) a étudié un modèle ana¬ 
logue en subdivisant les populations en plusieurs groupes d'âge. Les 
résultats obtenus sont plus convaincants mais restent soumis aux mêmes 
critiques : des caractères variables sont considérés comme constants. 

5. Relations entre proies et prédateurs. 

Equations de Lotka-Volterra 

Dans ce modèle, les relations entre les populations sont d’un type 
différent : on considère des espèces prédatrices (ou parasites) et leurs 
proies. 

Si H représente la population des proies et P les prédateurs (ou les 
parasites), on peut supposer que (Volterra, 1926, 1937, 1938; Lotka, 
1907. 1912, 1920, 1923, 1932 a et b) : 

dH 

-= (a, - b, P) H 

dt 

dP 

-= ( — a, + b. H) P 

dt 

a, = coefficient de reproduction de l'espèce 1 (proies); 

b , = coefficient de mortalité de l’espèce 1 due à l'attaque des pré¬ 
dateurs; 

a. . = coefficient de reproduction de l'espèce 2 (négatif = de mortalité, 
quand les proies sont absentes); 

b. . = coefficient de reproduction de l’espèce 2 quand des proies sont 
disponibles. 

On a : 

c/H (a, — b, P) H 

dP ~ (-a, + b.. H ) P 

dH dP 

a.. (-) — b.. dH + a, (-) — b , c/P = O 

H P 


D'où : 
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En intégrant, on obtient : 

a : lnH-i 2 H + a 1 lnP-6 1 P = C 

Cette équation peut se représenter par une courbe fermée, les varia¬ 
tions de H et de P étant cycliques. Pour des cycles de faible amplitude, 
la courbe résultante tend vers une ellipse (Chiang. 1954). 

Une équation aussi simple, donnant des cycles à allure si « biologique » 
a eu un très grand succès. Elle est très facile à simuler sur un ordinateur 
analogique (Anon., 1965). Thompson (1922 a. b, c. 1924, 1928, 1929) 
a largement utilisé ce type d'équation pour l'étude des parasites entomo- 
phages. Bodenheimer (1938) et Bodenheimer et al. (1952) ont pu 
mettre en évidence des variations des cycles en imaginant des actions 
catastrophiques. Mais ces manipulations (Utz, 1961; Utz et al.. 1963) 
ne doivent pas faire oublier que ces modèles sont basés sur des notions 
très discutables : 

— on imagine que les coefficients de reproduction et de mortalité sont 
indépendants de la densité de population. C'est évidemment une 
erreur, largement démontrée (Holling, 1959 a et b; 1961, 1963; Utida. 
1943; Watt. 1959); 

— on imagine aussi que — dans le cas du parasitisme — les parasites ne 
font pas la distinction entre les hôtes déjà parasités et les individus 
sains. Biologiquement, cette supposition est erronée : la plupart des 
femelles de parasites distinguent parfaitement les hôtes adéquats 
(Jones, 1937; Lloyd. 1938; Salt, 1932. 1934; Ullyett. 1936. 1943, 
1945, 1949 a et b). 

Des améliorations ont été apportées aux équations originales de Lotka- 
Volterra (Bailey, 1931. 1933 a. b. 1934; Nicholson et al., 1935: 
Odum et al., 1954; Varley, 1959; Wangesky et al., 1956. 1957) en fai¬ 
sant intervenir notamment des décalages dans le temps entre le moment 
de la fécondation et l'apparition des nouvelles générations. Ces additions 
restent malheureusement toujours très en deçà de la réalité biologique. 
Aucune de ces améliorations ne faisant intervenir par exemple que l'effi¬ 
cacité des parasites diminue lorsque la population augmente, ni l'influence 
des facteurs physiques pourtant bien connus (Burnett, 1953 a. b, 1956. 
1958. 1959; De Bach et al., 1941, 1945. 1947; Miller, 1959. 1960; 
Rescigno et al., 1965. 1967; Rescigno. 1968; Ullyett. 1949 a et b: 
Watt, 1959). 

6. Modèle à multi-relations. Méthode de Garfinkel et al. 

Les équations étudiées précédemment sont basées sur un petit nombre 
de relations entre les organismes. Depuis quelques années, l'utilisation 
des ordinateurs permet d'utiliser des modèles constitués d'un grand 
nombre d'équations simples (Chance et al., 1960, 1962; Garfinkel et al.. 
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1961, 1964 a, b, 1967; Garfinkel. 1962 a. b, 1965; Larson et ai, 1962; 
voir aussi Mann, 1969). 

La base des relations est une analogie avec la loi d’action des masses : 

A -f- B — * C + D 

- dA — dB 

-=- = k{ A) (B) 

dt dt 

Le modèle biologique considère l'existence des catégories biologiques 
suivantes : 

GRS = herbe; 

RTS = racines, inaccessibles aux herbivores; 

FRB = herbivore nourri, pouvant se reproduire; 

SRB = herbivore sous-alimenté, ne se reproduit pas. 

Des équations différentielles calculées numériquement sont associées : 

d (GRS) 

-= k, (RTS) + k, (GRS) - h (GRS)-' - k : (FRB) (GRS) 

dt 

- k, (SRB) (GRS) 
d (SRB) 

-= k< (FRB) - k s (SRB) - k- 2 (GRS) (SRB) 

dt 

d (FRB) 

-= k 3 (FRB) + k 2 (GRS) (SRB) - k 0 (FRB) - k, (FRB) 

dt 

Ces équations, qui peuvent être complétées par l'addition de prédateurs, 
demandent un certain nombre de remarques. 

Comme dans les modèles précédents, il est admis implicitement que 
les animaux se reproduisent constamment et de façon linéaire, et que 
les « jeunes » sont immédiatement aptes à se reproduire, ce qui est faux 
évidemment. De plus, les coefficients de reproduction sont supposés être 
indépendants de la densité de la population, ce qui est aussi inexact. 
Néanmoins, cette méthode est intéressante car elle permet d’introduire 
un grand nombre de paramètres. 


II. LES MODELES STOCHASTIQUES 

Les modèles déterministes décrits précédemment peuvent être modifiés 
et transformés en modèles stochastiques (Bailey, 1950, 1957, 1968; 
Bartlett, 1955, 1956, 1957, 1960; Bartlett et al., 1960; Bertrand, 
1962). 
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Diverses équations ont été modifiées. Nous n’en retiendrons que les 
principales. 

1. Population croissant dans un espace limité. 

Equation modifiée de Verhulst 

Les coefficients de reproduction et de mortalité du modèle détermi¬ 
niste sont remplacées par : 

la probabilité d'une naissance au temps dt est (<*,N — /?,N-) dt; 

la probabilité d'un décès au temps t est (a_,N — /LN 2 ) dt. 

Les coefficients sont réunis par la relation : 

ai — a,. = a 

fil + & = P 

Si dt est petit, on peut écrire (Bartlett, 1957) : 

dN = (aN — /3N-) dt + dZ, — dZ 2 (1.1) 

dZ, et dZ z étant des variables indépendantes de Poisson à moyenne 
nulle et de variance : 

pour dZ, : (a,N — /?,N 2 ) 

pour dZ_, : (a 2 N — /? 2 N 2 ) 

La comparaison entre les résultats obtenus en utilisant les modèles 
stochastiques et déterministes de ce type a été effectuée (Feller, 1939. 
1940; Bailey, 1950). Les valeurs calculées avec l'équation (1.1 ) semblent 
se rapprocher davantage de la réalité. Mais elles restent d'un usage 
très difficile. Dans la pratique, on ne voit pas très bien ce que ce genre 
d’équations apporte de nouveau quant à la compréhension des phéno¬ 
mènes biologiques. 

2. Relations entre proies et prédateurs. 

Equations modifiées de Lotka-Volterra 

L'équation de Lotka-Volterra : 
dN, 

-= (a, — b u — b r . N 2 ) N, 

dt 

peut s'écrire de façon plus précise (PiÉLOU, 1969) : 

dN, 

-= N. (A, (N,/N 2 ) - (N,/N 2 )) 

dt 
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où A, (N,/N-J est le coefficient de reproduction de N, en présence de Nj 
représentants de sa propre espèce et de N, représentants de l'espèce 2. 
et où /i, (N, /NJ est le coefficient de mortalité de N, en présence de N, 
représentants de l'espèce 1 et de N., représentants de l'espèce 2. 

Si 

A, (N,/NJ = a, b - 6,, b N, - 6 a2b N 
Mi (N,/NJ = a, d + f>u d N, + 6 Ud N, + f> l2(1 N., 

On a : 

c/N, 

—-=N. ((«.„-a ld ) - (/>,.„-/>„„) N, - NJ 

at 

c/N, 

— = N, ((a 2 „-a 2 J - (f» 21b - fc sld )N, - (/>„„-6 22(I )NJ 

dt 

En utilisant ces équations et des valeurs de coefficients variables, 
Barnett (1962) a établi un certain nombre de courbes. Mais encore une 
fois, la méthode pèche par la base : les coefficients choisis n’ont pas 
de réalité biologique stricte. Les faire varier au hasard ne fait que masquer 
leur complexité sans les résoudre en éléments plus simples. 

D'autres modifications aux équations de Lotka-Volterra ont été 
apportées par divers auteurs en substituant des coefficients stochastiques 
aux coefficients déterministes mais sans résultats nouveaux appréciables 
(Bartlett, 1949, 1957, 1960; Chiang. 1954. 1957; Kendall. 1948, 1949; 
Leslie, 1957, 1958. 1959, 1962; Leslie et al., 1958). 

En conclusion, la plupart des modèles que nous venons d’étudier pré¬ 
sentent certains défauts. Leurs hypothèses de base sont trop générales 
et envisagent des situations biologiques complètement inexistantes. 

Le modèle que nous proposons tente d'éviter ces critiques en simulant, 
de la façon la plus simple possible, un certain nombre de relations bio¬ 
logiques de base. 


MODELE MATHEMATIQUE PROPOSE 
I. METHODES ET PARAMETRES 
1. Choix des paramètres 

La validité d'un modèle théorique dépend évidemment du choix des 
paramètres utilisés. Pour être réaliste, un modèle doit posséder un 
grand nombre de variables indépendantes et doit être structuré de façon 
telle que les facteurs impliqués interviennent au bon moment, aussi bien 
dans le temps que dans l’espace (Watt. 1955, 1956. 1959 a, b, 1960, 
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1961 a. b, c, 1962, 1963, 1966. 1968). Le modèle doit aussi décrire 
les relations complexes qui existent entre les différents paramètres bio¬ 
logiques et physiques. Même si certains facteurs n'ont pas encore été 
suffisamment étudiés sur le terrain, par exemple l'effet des variations 
climatiques sur la répartition et la survie d’une espèce, ils doivent être 
prévus dans les équations. Ces manques ne sont, en effet, que temporaires 
et sont quantifiés au fur et à mesure des recherches (Holling, 1959, 
1961; Mook et al.. 1960; Prop. 1960: Tinbergen. 1960; Dixon, 1959; 
Kramer, 1960). 

La plupart des modèles proposés ne font pas intervenir les variations 
des conditions physiques ou chimiques (Beverton et al., 1957; Nicholson 
et al., 1935; Silliman, 1966; Walters. 1969) ce qui amoindrit considé¬ 
rablement leur intérêt. 


a. Paramètres physiques 

Dans le modèle présenté ici. différentes variations climatiques ont été 
incluses. On ne considérera cependant, pour commencer, que 2 variables 
fondamentales : la température et l'humidité relative. Ces deux para¬ 
mètres peuvent prendre n'importe quelle valeur suivant les nécessités 
de la simulation. Le fait important est que ces conditions climatiques 
sont rendues indépendantes dans deux zones différentes. Au même 
moment, deux milieux physiquement différents sont disponibles pour les 
organismes. Ce sont leurs paramètres biologiques, et eux seuls, qui vont 
être responsables de la répartition (et donc, en partie, de la densité). 

b. Paramètres biologiques 

Un milieu écologique, même très simple, est peuplé de différentes 
espèces vivantes ayant chacunes des caractéristiques différentes variant 
en fonction du temps, de leur état physiologique (reproduction, etc.) 
et des conditions climatiques. 

Fondamentalement, on peut distinguer des organismes vivants immo¬ 
biles (les végétaux et certains stades de la vie animales : œufs, larves, 
animaux endormis) et des organismes mobiles (animaux). Les relations 
entre les organismes d’un milieu écologique font donc intervenir une 
notion de déplacement caractérisé par une vitesse et un sens. 
Envisageons cette première notion. 

— Caractéristiques du déplacement des êtres 
vivants 

Ce paramètre — déplacement des organismes — est évidemment fon¬ 
damental dans la mise au point d'un modèle. Du très grand nombre des 
observations biologiques existantes, il faut déduire des relations générales. 
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Comment, par exemple, un animal trouve-t-il sa nourriture ? 

Certains parasites de Diptères (Alysia mandacutor Panz. et Mormo- 
miella vitripennis Walk.) sont attirés par les cadavres. En effet, leurs 
proies se trouvent la plupart du temps occupées à dévorer les charognes 
(Laing, 1937). Le principe attractif est l'odeur de putréfaction et notam¬ 
ment les molécules d'indol et de skatol (Hobson. 1935. 1936). 

Les femelles de Drosophila pondent leurs œufs sur les fruits fer¬ 
mentés et sont attirées par l'odeur de fermentation (Barrows. 1907). 

Cette réceptivité aux odeurs peut varier avec l'état physiologique 
(Thorpe ef al., 1939). 

La vue est soumise à des variations semblables. Le Papillon Pieris 
brassicae L. est attiré par la couleur des fleurs quand il est affamé mais, 
en période de ponte, il est spécialement réceptif à la couleur verte des 
feuilles (Ilse. 1937). 

D'autres animaux combinent plusieurs sens dans une même approche 
(Tinbergen, 1932. 1935; Ullyett et ai, 1940; Warnke, 1931). 

Comment arriver à quantifier ces différents comportements ? 

Divers auteurs (Varley et al., 1963. 1968; Hassell, 1969) ont ainsi 
défini, dans le cas des parasites, une surface de découverte (a) exprimant 
la possibilité pour une femelle de découvrir et d'attaquer un hôte adéquat. 
Si la densité de la population de parasite est p, la densité initiale des 
hôtes est h„ et la densité finale des hôtes est h„ on a : 

K 

p . a = log e - (1.1) (Hassel, 1969) 

hr 

Nicholson (1933) et Nicholson et al. (1935) ont supposé que cette 
surface de découverte diminue lorsque la densité des parasites augmente. 

Mais d'autres facteurs peuvent influencer (a) : 

Chez Caccellia gnava (Meig.) (Tachinide). la recherche des larves 
Bupalus pinarius (L.) s'intensifie avec l'ensoleillement (Klomp, 1959). 

Chez l'ichneumon Nemeritis canescens (Grav.), à la recherche des 
chenilles d'Ephestia cautella (Walk.), la surface de recherche augmente 
lorsque la densité des hôtes augmente (Takahashi, 1968), atteint un 
maximum puis se réduit (Holling, 1959). 

A la recherche des Anagasta kühniella (Zell.), en laboratoire, la 
valeur de (a) peut se quantifier suivant : 

log, 0 a = — 1.22 — 0.54 1og,„p 

(Huffaker et al., 1969) 
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En reprenant les résultats de différents auteurs (Burnett, 1956, 1958; 
Ullyett, 1949 a, b), Hassel et al. (1969) ont pu définir (a) suivant : 

log 10 a = log 10 Q — m log„, p 

avec Q = coefficient de recherche; 

m = coefficient d’interférence mutuelle. 

En remplaçant dans (1.1), on obtient : 


Q 



h 0 pQ 

log,.-=-= p*' m Q 

h, p» 1 

Hassell et al. (1969) ont montré que plusieurs espèces de parasites 
possédant des valeurs différentes de m (coefficient d'interférence 
mutuelle) et Q (coefficient de recherche) peuvent parfaitement coexister 
dans le même milieu. 

Ce paramètre (a) englobe, en réalité, deux notions biologiques : une 
vitesse de déplacement et un diamètre de perception. Dans le modèle que 
nous présentons, vitesse et diamètre sont considérés séparément et varient 
suivant les espèces, leur état physiologique et leur stade de développe¬ 
ment. Pour l’instant, nous ne considérons pas de facteur m; il sera intro¬ 
duit dans la suite. 

Tous les modèles proposés jusqu'ici suppose que le déplacement des 
animaux s'effectue au hasard. Dans le cas des Mammifères, cette hypo¬ 
thèse n’est pas exacte (Siniff et al., 1969) notamment par la présence 
d'un territoire préférentiel (Calhoun et al., 1958; Stumpf et ai, 1962); 
cette surface de territoire peut être quantifiée (Cagle, 1953, 1956; Fitch. 
1949; Jennrich et al., 1969; Skutch, 1950). Une circulation au hasard 
est cependant vraisemblable dans le cas des Invertébrés. Dans le modèle 
que nous proposons, le déplacement des organismes est supposé s'effec¬ 
tuer au hasard mais le vecteur exprimant le déplacement n'est pas linéaire 
mais assimilé à une surface équivalente à un cercle dont le rayon est 
défini pour chaque espèce. Dans le cas d'un prédateur, le rayon exprime 
la sensibilité de ses organes de perception; dans le cas d'une proie, le 
rayon exprime son influence sur l'environnement (odeur, bruit, etc.). 
Quand deux de ces surfaces sont jointives, il est évident que le déplace¬ 
ment du prédateur ne s'effectue plus au hasard : il se dirige immédiate¬ 
ment vers sa proie. Le modèle concilie donc les deux aspects de base du 
comportement : recherche avec une certaine vitesse et certaines possi¬ 
bilités de perception. 

L’équation utilisée dans ce cas a été définie et discutée précédemment 
(Tercafs, 1968). Elle exprime le nombre de rencontres entre les indi¬ 
vidus de n populations en fonction des paramètres définis ci-dessus. 
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Dans le cas des individus jeunes, elle est tempérée par un paramètre 
supplémentaire e qui tient compte du fait que les jeunes prédateurs ne 
mangent pas toutes les proies qu'ils rencontrent : ils effectuent un choix 
(ce paramètre e est équivalent au paramètre w i} de O’Neill (1969)). 

— Caractéristiques de la reproduction 

Les caractéristiques de la reproduction que nous avons retenues sont 
assez élémentaires mais la plupart des modèles mathématiques ne les 
envisagent même pas (voir cependant Leyton, 1968; Talus et al., 1969 
et Taylor, 1968 mais où les naissances ne s'effectuent pas de façon déter¬ 
ministe). 

Nous considérons que les animaux sont constitués de mâles et de 
femelles dans une proportion définie. Pendant une certaine période, la 
rencontre entre un mâle et une femelle de même espèce provoque un 
accouplement. La femelle, qui ne peut être fécondée qu'une seule fois, 
subit alors une période de gestation menant à la ponte des œufs. Ces 
œufs sont groupés en ponte qui peuvent être plus ou moins protégées. 
Le développement des œufs se poursuit et ils donnent naissance à des 
formes larvaires dont les caractéristiques biologiques vont varier en fonc¬ 
tion du temps. Par rapport aux adultes, leur vitesse de déplacement, 
leur diamètre de perception ou d'influence, leur possibilité de dévorer 
des proies, etc, sont différentes. 

Le modèle ainsi défini rencontre donc un très grand nombre de carac¬ 
téristiques biologiques de base habituellement regroupées (Legay, 1968). 

— Caractéristiques de la répartition 

Des milieux proches ne possèdent pas nécessairement des conditions 
climatiques identiques. Suivant leurs préférendums, les animaux peuvent 
se répartir différemment dans diverses régions. Les variations clima¬ 
tiques vont entraîner des modifications dans la densité des populations 
d'un milieu. Ce fait est évidemment ignoré dans un modèle simple où 
l'on considère un climat stable. Dans le modèle que nous présentons, 
nous avons ajouté des paramètres climatiques et corrélativement des 
préférendums aux animaux. Dans un premier stade, nous allons consi¬ 
dérer uniquement la température et l'humidité relative. 

— Caractéristiques de la nutrition 

Les chaînes trophiques, même les plus complexes, ne sont en fait 
constitués que d’une relation simple : un être vivant ou un groupe sert 
de proie à un autre groupe (Mac Arthur, 1955; Reddingius, 1963). 
Cette relation est fondamentale et toutes les autres en découlent. Pour 
ne pas compliquer exagérément le modèle, nous avons ainsi choisi deux 
groupes de « carnivores » et 4 groupes « d'herbivores ». Les premiers 
se nourrissent au dépens de seconds. Leur capacité de destruction est 
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limitée à un certain nombre de proies par unité de temps. Les carnivores 
qui ne sont pas nourris sont dits « sous-alimentés » et des caractéristiques 
physiologiques différentes ( = des paramètres de valeurs différentes) 
leur sont affectés. Les herbivores, par contre, possèdent des sources de 
nourriture en quantité nécessaire, en permanence. Ultérieurement, il 
suffira de compléter le modèle en définissant une norme de croissance 
et une répartition des végétaux. Les stades larvaires de chaque catégorie 
possèdent des caractéristiques différentes de celles des adultes. 

Ces paramètres biologiques étant ainsi définis, nous allons établir, 
entre eux, un certain nombre de relations. 

2. Relations entre les paramètres 
a. Pref etendum 
— Thermopreferendum 

D'une façon générale, le thermopreferendum d'un groupe animal, dans 
un état physiologique donné et à un moment donné, peut être représenté 
par une courbe. 

Dans un cas idéal où, toute autre chose étant égale, une vaste gamme 
de température est présentée, les animaux vont se distribuer suivant leur 
température préférentielle. Cette distribution peut prendre différentes 
formes, mais avec trois caractéristiques fondamentales. Chaque preferen- 
dum est caractérisé par une température minimale (ANLOT) (1) et par 
une température maximale (ANLTT). En dehors de ces extrêmes, les 
conditions ne sont plus favorables à la vie de l'animal. En outre, on peut 
distinguer une valeur maximum où la plupart des individus du groupe se 
rassemblent. Ce maximum est donné par la valeur de ANVLMT. 

Dans le modèle défini plus loin, le thermopreferendum sera représenté 
par une courbe sinusoïdale et l'intervalle séparant les températures préfé¬ 
rentielles extrêmes sera divisé arbitrairement en 100 parties égales. Cette 
disposition n’est choisie qu'à titre d’exemple et peut être modifiée suivant 
les impératifs des travaux. 

Dans le cas choisi ci-dessus, les coefficients de répartition R T ° en 
fonction de la température seront donc de la forme : 

R t » = ANVLMT X sin « (1.1) 

(avec a variant de O à tt) 

Dans les calculs définis par les programmes FORTRAN décrits plus 
loin, la formulation est plus complexe car il faut au préalable calculer 
100 valeurs particulières de l’abscisse et y associer les 100 valeurs définies 
par l'équation (1.1). 

( 1 ) Les abréviations des paramètres sont identiques à celles utilisées dans les pro¬ 
grammes FORTRAN et définis page 47. 
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— Hygroferendum 

D'une façon analogue au thermopreferendum, on peut définir des 
coefficients de répartition en fonction de l’humidité relative. 

La formulation est semblable à celle définie par l'équation (1.1) et est 
de la forme : 

R r „ = ANVLMH X sin « (2.1 ) 

(avec a variant de 0 à tt) 

— Conflits entre preferendums. Les valeurs 
g a m m a s 

Les notions de preferendums définies ci-dessus et qui sont les défini¬ 
tions classiques, n'envisagent en aucune façon le cas de conflits entre 
preferendums. Quelle est l’origine possible de ces conflits ? Un exemple 
fera mieux comprendre. Soit un groupe animal AN possédant le thermo- 
préférendum simplifié : 

R t „ = ANVLMT = AN (3.1) 

(avec T° = Tl ) 

Ce qui revient à dire que tous les animaux du groupe se localisent 
dans la zone possédant la température Tl. 

Supposons, d'autre part, que ces animaux AN possèdent un hygro- 
préférendum de la forme : 

Rk„ = ANVLMH = AN (3.2) 

(avec RH = RH1) 

Ce qui veut dire que tous les animaux AN se groupent, quand ils en 
ont le choix, dans la zone possédant une humidité relative RHl. 

Plaçons maintenant ces animaux dans une enceinte expérimentale possé¬ 
dant deux zones communicantes possédant les caractéristiques physiques 
suivantes : 

Zone 1 : température : Tl 

humidité relative : RH2 

Zone 2 : température : T2 

humidité relative : RHl 

T2 et RH2 étant très différents de Tl et RHL 

Les animaux AN qui dans l'hypothèse de départ recherchent un milieu 
possédant les conditions physiques Tl et RHl se voient donc obliger 
de choisir entre les deux milieux présentés. Ce choix, obligatoire puis¬ 
qu'il y a conflit entre les preferendums, peut être quantifié si on introduit 
une notion nouvelle, celle d'importance relative accordée par 
l'animal à chacun de ses preferendums. Nous appelerons cette notion : 
gamma (y). 
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Dans l’exemple précédent, si les animaux accordent une importance 
relative très élevée au thermopréférendum et négligent l'hygropréfé- 
rendum, ils vont tous se retrouver dans la zone 1. Mais un cas aussi 
absolu ne se présente jamais; c’est pourquoi il faut distinguer deux 
valeurs de gammas : 

— un gamma thermique (y T .) exprimant l’importance relative accordée 
par les animaux à leur thermopréférendum; 

— un gamma hygrométrique (y RH ) exprimant l'importance relative accor¬ 
dée par les animaux à leur hygropréférendum. 

Ces gammas peuvent s'exprimer en %. La répartition des animaux AN 
considérés plus haut est donc sous la dépendance de ces divers facteurs 
suivant l'équation (3.3) : 

Rp = (yT X R-i-) + (Yrii X Rrh) (3.3) 

avec R ( , = répartition totale. 

Un exemple permettra de visualiser ces importantes notions : 

Considérons un groupe animal AN. possédant les caractéristiques sui¬ 
vantes : 

— thermopréférendum : 

AN LOT = 5. 

ANLTT = 20. 

yi = 80. 

— hygropréférendum : 

ANLOT = 40. 

ANLTH = 100. 

y R H = 20. 

Rkh <0 . m = 25 % 

Rk.w„ = 50 % 

Rh„ s0 . 100 = 25 % 

Les animaux ont le choix entre deux zones expérimentales (ZE et ZP) 
de caractéristiques différentes : 

ZET = 12 ZPT = 13 

ZERH = 82 ZPRH = 85 


Compte tenu de leurs caractéristiques, les animaux AN vont se répartir 
entre les deux zones. On a. suivant l’équation (3.1) et (3.2) : 


Dans la zone ZE : 


Rt» 

Rrh 


= 50 
= 25 


R T-. 

Rrh 


= 50 
= 25 


Dans la zone ZP : 
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La répartition générale est, suivant l'équation (3.3) : 

Pour la zone R,, = 90 
ZE 

Pour la zone R,, = 90 
ZP 

Les deux coefficients de répartition étant identiques les animaux se 
distribuent à raison de 50 % dans chacune des deux zones. 

Mais si les caractéristiques physiques des zones ZE et ZP varient, 
la répartition va être tout à fait différente. En effet, soient les conditions 
suivantes : 

ZET = 8 

ZERH = 75 

Les coefficients de répartition deviennent : 

Pour la zone ZE : R p = 50 

Pour la zone ZP : R p = 90. 

Ramenés en % suivant l'équation (3.4). on obtient : 

ZE R, = 35.7 
ZP R„ = 64.3 

Des modifications dans les conditions physiques provoquent donc des 
remaniements dans la répartition des animaux. 

Cette notion de gamma doit cependant être modifiée pour se rappro¬ 
cher davantage de la réalité. En effet, l'importance relative qu’un animal 
accorde à un préférendum varie avec les circonstances. Ce qui nous 
amène à considérer des gammas variables. 

D'une façon générale, on peut considérer que l'animal va accorder une 
importance de plus en plus grande à un préférendum plus on s’écarte 
des conditions optimales. 

Dans le modèle défini plus loin, nous avons considéré, de ce fait, que 
les gammas (y T et y R n) pouvaient se représenter par les équations sui¬ 
vantes : 

y T = ANRLMX + (ANRLMI X cos «) avec a variant de ir à 3w 

2 2 

(3.5) 

(3.6) 


y R h se divise en deux parties : 

'/ R» (i AXRHO) “ ANRELM 


ZPT = 14 
ZPRH = 100 


y ru 


(AXRHO + 1 ■ 


Iu0) = (ANRELM - ANREMI) o/k-axrho. + ANREM1 

(3.6) 

avec K variant de ANRHO + 1 à 100. 
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A chaque température et à chaque humidité relative est donc associée 
une valeur particulière de gammas. 

L'ensemble de ces notions de préférendums et de gammas va être 
utilisée dans le calcul de la répartition des animaux dans un milieu 
à n zones écologiquements différentes. 

b. Répartition des animaux 
dans un système à n zones relativement différentes 

— Généralités 

Considérons un système écologique constitué de différentes zones join¬ 
tives possédant des conditions écologiques (humidité, température) parti¬ 
culières : Tl, RH1, T2. RH2, T3. RH3. 

Des animaux AN sont répartis dans le milieu 1. Imaginons que leurs 
organes des sens leur permettent de percevoir des variations spatiales des 
conditions physiques. Si les animaux AN parcourant le territoire 1, ils 
vont rencontrer, à un certain moment, la zone de contact entre les milieux 
1 et 2. Dans cette zone, les variations locales des conditions physiques 
donnent à l'animal la possibilité d'un choix : ou bien il préfère les condi¬ 
tions de la zone 1 ou bien il préfère les conditions de la zone 2. Dans ce 
deuxième cas, il entre dans la zone 2 et commence à la parcourir et, 
comme précédemment, peut entrer en contact avec la zone 3. Le même 
processus peut se répéter ainsi indéfiniment. Mais, comme on le voit, 
les animaux n’ont, à chaque fois, qu’à effectuer un choix qu'entre deux 
possibilités. Dans un modèle, il suffira donc de mettre au point la méthode 
permettant de calculer la répartition des animaux entre deux zones uni¬ 
quement. et de répéter les calculs suivant le nombre de zones envisagées. 
Dans le modèle défini plus loin, nous n'envisagerons que deux zones, pour 
simplifier cette première approche. 

— Méthodes 

Rappelons d’abord que dans les préférendums décrits précédemment 
définissent une distribution idéale que l'on peut subdiviser en intervalles 
variables suivant le degré de perception des animaux. Le nombre d'inter¬ 
valles doit être plus élevé si les animaux sont capables de distinguer les 
variations minimes dans les conditions physiques. 

Dans la zone de température où le thermopréférendum existe, on 
définit donc un certain nombre de coefficient de répartition. Le nombre 
est fixé à 100 dans le modèle proposé. On a : 

ANLTT 

V RPANT = 1 (4.1) 


A N LOT 
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La somme des coefficients de répartition-température des animaux AN 
compris dans l'intervalle où le thermopréférendum existe (intervalle 
ANLOT-ANLTT) est égal à 1, par définition. 

De même, dans le cas de l'hygropréférendum, on a : 

ANLTT 

V RPANH = 1 (4.2) 

A N LOT 

Le milieu étudié est constitué de deux zones écologiquement différentes. 
Dans la zone 1, la température est de ZET "C et l'humidité relative est 
de ZERH %RH. Dans la zone 2, la température est de ZPT "C et 
l'humidité relative est de ZPRH % RH. 

Plusieurs cas peuvent se présenter suivant les valeurs relatives des 
conditions physiques. Considérons d'abord le cas de la température. 


1 er cas : les 2 zones ont des conditions identiques. 

ZET = ZPT 

Dans ce cas extrêmement simple, les 2 zones n’en forment plus qu’une 
au point de vue température et on est ramené au cas simple précédem¬ 
ment défini (4.1 ). 


2"" f cas : Les 2 zones ont des conditions phy¬ 
siques différentes (cas général) 

ZET ^ ZPT 

Considérons d'abord la température moyenne résultante. On a : 
TMOY = (ZET + ZPTJ/2. 

Le problème qui se pose est de savoir ce que font les animaux quand 
ils ont le choix entre des températures également distantes de la tempé¬ 
rature optimale. Nous allons considérer que les animaux se distribuent 
de la façon suivante : 

— une partie des animaux se rendent dans la zone 1 suivant : 

TMOY 

RTE = V RPANT 


A N LOT 
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— l'autre partie se rend dans la zone 2 suivant : 

ANLTT 

RTP = V RPANT 
TMOY 


RTE et RTP expriment donc la répartition des animaux considérés 
dans les 2 zones, au moment étudié. 

Un raisonnement analogue s'applique dans le cas de l'humidité rela¬ 


tive. On a : 


dans la zone 1 


H MO Y 

RHE = V RPANH 
ANLOH 

ANLTH 


dans la zone 2 RHP = ^ RPANH 

HMOY 

HMOY étant l'humidité relative moyenne. 


Ces particularités fondamentales de la répartition étant définies, il nous 
faut maintenant faire intervenir les coefficients gammas définis précédem¬ 
ment. A un moment défini Tl, deux gammas sont définis : 

^T(ti) = A 


Y hr,td = B 

Les valeurs A et B sont données par les caractéristiques respectives 
des 2 gammas. Mais ce qui importe dans le calcul ultérieur de la répar¬ 
tition, ce ne sont pas les valeurs des gammas en tant que telles mais leur 
importance relative. On calcule des gammas globalisés: 

— gamma globalisé pour la température : 

I'tcmti) = ( t 11 / ( Y t i ti ) + Vriiitd) 

— gamma globalisé pour l’humidité relative : 

*7rhg(ti) = '7rh(ti)/(’>'t(ti) + Trimtd) 

Chacun des gammas ainsi calculé est associé au coefficient de répar¬ 
tition correspondant. On obtient : 

Pour la zone 1 : 

ZERP (Tl) = (RHE x Ytg,t, >) + (RHE X riigitd) 
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Pour la zone 2 : 

ZPRP(Tl) = (RTP X 7 to(td) 4* (RHP X Ykiigiti)) 

ZERP et ZPRP étant les coefficients de répartition dans chacune des 
2 zones. Les valeurs finales sont alors calculées proportionnellement : 

ZERP(Tl) = ZERP(T1 ) / RP 
ZPRP(Tl) = ZPRP(T1 ) / RP 
avec RP = ZERP(Tl) + ZPRP(Tl). 

Si des animaux AN ont le choix entre 2 zones écologiquement diffé¬ 
rentes (1 et 2), leur répartition, au temps Tl, sera : 

- dans la zone 1 : AN (1) = AN X ZERP(Tl) 

— dans la zone 2 : AN (2) = AN X ZPRP(Tl). 

C. Activité 

La première partie exposant les relations entre les paramètres a défini 
les facteurs régissant la distribution des animaux dans des milieux écolo¬ 
giquement différents et ce, en considérant les caractéristiques les plus 
importantes de leur comportement vis-à-vis de ces facteurs. 

Il nous faut maintenant étudier les autres aspects du comportement 
intervenant dans 1 activité des animaux. Comme signalé précédemment 
(page 17), nous allons définir dans notre modèle 2 grandes catégories 
d'animaux : les carnivores et les herbivores. Il est bien évident qu'il ne 
peut y avoir interactions entre ces deux groupes d'animaux que lorsqu'ils 
sont en contact. Ce nombre de rencontres peut être quantifié. 

L'équation que nous avons utilisée a été définie dans un travail anté¬ 
rieur (Tercafs, 1968) : 

m 1 Ç 2 -r f - 

0, 2 = n . — . (c/l + cf2) I | V\ + V 2 2 - 2V,V- cos« d« 

U 2,j, 

où : 

0,., = nombre de rencontres entre les populations 1 et 2; 
n = nombre d'individus dans la population 1; 

m = nombre d'individus dans la population 2; 

U = surface du milieu; 

dl = diamètre de perception ou d'influence d'un individu 1; 

d2 = diamètre de perception ou d'influence d'un individu 2; 

V1 = vitesse moyenne des individus 1 ; 

V2 = vitesse moyenne des individus 2. 
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Dans le modèle mathématique proposé, il y 4 catégories de carnivores 
et 4 catégories d'herbivores. Nous nous proposons d’étudier l'évolution 
de ces différentes populations en fonction du temps. Ce temps sera défini 
par un indice (K). 

Les caractéristiques des populations étudiées sont définies de la façon 
suivante (tableau 4.1). 

TABLEAU 4.1. 


Caractéristiques des populations étudiées. Symboles utilisés 


Catégorie 

Etat 

Nombre 

Diamètre 

Vitesse 

carnivore . 

nourri 

AN 

DA 

VA 

carnivore . 

nourri 

BN 

DB 

VB 

carnivore . 

sous-alimenté 

AS 

DA 

VAS 

carnivore . 

sous-alimenté 

BS 

DB 

VBS 

herbivore . 

nourri 

C 

DC 

VC 

herbivore . 

nourri 

E 

DE 

VE 

herbivore . 

nourri 

F 

DF 

VF 

herbivore . 

nourri 

G 

DG 

VG 


Le nombre de carnivores nourris pendant l'unité de temps est défini 
par le nombre de rencontres entre ce groupe de carnivore et les groupes 
d'herbivores vivant dans le même milieu. Si l'on considère la population 
AN, par exemple, on aura : 

C ( K ) 1 * 2"" 

0wo = AN(K) .-.(DA 4- DC) .— i *’ Y VA* 4- VC 2 - 2VA VC cos «da 

U 2 tt J 0 


9 AN K 


E(K) 1 

AN(K) .-. (DA 4- DE) .— 


U 


2 7T 


. <£7T , - 

\ Y VA 2 4- VE- - 2VA VE cos« de 


F(K) 1 r 2r - 

tf WF = AN(K).-.(DA + DF).— \ VA 2 4- VF- — 2VA VF cos« da 

U 2tt J 0 


G(K) 1 

*ano = AN(K) .-. (DA + DG) .— 

U 2 TT 


r 2 - ,- 

^ Y VA 2 4- VG 2 - 2VA VG cos« da 


avec : 


AN (K 4- 1) = ^ ANC 4- 0 ANE 4- 0ANF 4" 9AÜO 


(5.1) 
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AN (K) défini ici le nombre de carnivores nourris pendant l'intervalle 
de temps considéré. Mais le nombre de carnivores sous-alimentés ( = qui 
n'ont pas mangé une proie pendant cette période) est donné par : 

AS (K) = AN (K - 1) - AN(K) + AS(K - 1) (5.2) 


De la même façon, on peut calculer le nombre d'herbivores restant 
dans le milieu après l'attaque des carnivores. On a : 


C(K) — C(K — 1 ) — 0AXCIK— 1 ) — 0ASCIK—1) — 0BNCIK— 1 ) — ^BSC(K-l) 


(5.3) 

0A8C. #bsc. ^iisc étant le nombre de rencontre entre les animaux AS, BS 
et BN avec les herbivores C. 

On déduit des équations analogues pour les autres catégories d'herbi¬ 
vores. 

Les équations définies ci-dessus ne sont valables que si le nombre de 
rencontres d'un carnivore avec les différentes catégories d'herbivores 
ne dépasse pas le nombre de carnivores présents dans le milieu. En d'au¬ 
tres termes, les équations ne sont applicables que lorsque les carnivores 
ne mangent qu'une proie par unité de temps. Envisageons maintenant le 
cas général. 


Si l'on considère l’équation 5.1, on voit que si AN(K) >AN(K—1 ) 
cela veut dire que le nombre de rencontre des carnivores AN avec les 
herbivores C. E. F et G a été plus grand que le nombre de carnivores 
présents. Les carnivores ont donc mangé plus d une proie par unité de 
temps. On a : 

PRMGAN(K) = AN(K)/AN(K — 1 ) 
avec PRMGAN(K) = nombre de proies mangées par unité de temps 
par les animaux AN. 


Tous les animaux AN étant nourris pendant cet intervalle de temps, 
il vient immédiatement : 

AN(K) = AN (K - 1) 


Les carnivores ne peuvent manger un nombre indéfini de proies par 
unité de temps; il existe une valeur maximum (PRGXAN) où leur 
appétit est saturé, c'est-à-dire au moment où PRMGAN(K) >PRGXAN. 

Les équations précédentes sont complétées de la façon suivante, où l'on 
définit des nouvelles valeurs 6' pour 6. 

SOMAN(K) = #ANC<K> + 0AXE(K> + ^ANP(K) + ^ANO(K) 

0 ANC(K) = (SOMAN(K)/0a N - C(K ,) .PRGXAN 

0'axe,k, = (SOMAN(K)/0ax E <k,) - PRGXAN 

0’anf,k, = (SOMAN(K)/<W,k.) .PRGXAN 

0’ang«k, = (SOMAN(K)/0 A n G ,k>) - PRGXAN 
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6 définissant le nombre réel de proies mangées et non plus 
le nombre de recontres. 

Accessoirement, les valeurs de chaque catégorie animale peuvent être 
modifiée par un facteur supplémentaire englobant des paramètres pas 
encore détaillés. Dans le modèle présenté, chaque catégorie est affectée 
d’un coefficient de mortalité exprimant les pertes subies pour des raisons 
autres que la prédation (maladie, vieillesse, etc.). 


d. Reproduction 


1. — Accouplement 

A un moment défini I débute la période d'accouplement. A ce moment, 
chaque population est considérée comme étant constituée de 2 sous- 
populations. les mâles et les femelles, la proportion de l'une et de l'autre 
pouvant être définie par un certain coefficient N. 

Pour chaque catégorie, on a donc : 

ANM (I) = AN(I) .NA 
ANF (I) = AN(I) . (1 —NA) 

ANM (I) = nombre de mâles de la catégorie AN au temps 1. 

ANF (1) = nombre de femelles de la catégorie AN au temps I. 


Le nombre d’accouplements est défini, à chaque unité de temps, par 
le nombre de rencontres entre les mâles et les femelles. On a : 


ANF(I) VA f 

ACCAN(I) = e AS manf = ANM(I) .- .DA.- \ 

U » J 


j/ 2(1 - 


cosrt) c/« 


Ce qui donne également le nombre de femelles fécondées = ACCAN( 1 ). 

Les femelles n'étant fécondables qu’une fois, il faut décompter les 
femelles déjà fécondées à chaque unité de temps : 

ANF (I + 1) = ANF(l) - ACCAN( 1 ) 

Les accouplements ont lieu pendant une certaine période (de 1 à 1, 
par exemple). A la fin de cette période, le nombre total de femelles 
fécondées (ACCANT) est donc de : 

II 

ACCANT = V ACCAN(I) 

I 
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On en déduit le nombre de femelles fécondées (ANFF) par rapport 
à la population totale : 

ANFF = ACCANT/ANF (II) 

— Ponte 

Après une certaine gestation, les femelles fécondées vont pondre des 
œufs. Le nombre moyen d’œufs pondus est caractéristique de l'espèce 
et de son état physiologique. Dans le modèle décrit, nous avons à nouveau 
subdivisé chaque population en 2 sous-populations de caractéristiques dis¬ 
tinctes. la proportion de l’une et de l’autre étant donnée par un facteur 
particulier (TRA). On obtient les 2 groupes suivants : 

i ZAN caractérisé par OAN 

AN 

r ZTRAN caractérisé par OTRAN 
avec ZTRAN = AN . TRA 
ZAN = AN . (1 - TRA) 

OAN et OTRAN caractérisant le nombre d'œufs pondus par une 
femelle de chaque catégorie. 

Les œufs sont pondus par groupe constituant une ponte, de diamètre 
DZA. Le nombre d'œufs dans chaque ponte est caractéristique du groupe 
et est donné par la valeur d'un paramètre (POAN). 

Les œufs appartenant à la catégorie ZTRAN sont affectés d'un coeffi¬ 
cient supplémentaire (PP) exprimant, que d’une façon ou d’une autre, 
ces œufs sont partiellement protégés de l'attaque des prédateurs. Ils 
peuvent être enterrés, camouflés, etc. 

Toutes les catégoéries d’œufs peuvent être détruites par les carnivores, 
pendant toute la durée de leur développement. L’expression de cette 
action est donnée par l'équation suivante : 

ZTRAN (K — 1 ) 

ZTRAN (K) = ZTRAN (K - 1) - ((- 

U 

. ( ( ( AN(K) + AS( K) ) . (DA + DZA) .VA) 

+ ( (BN(K) + BS(K)) . (DB + DZB) . VB))) . PP . BE) 

BE représente le nombre d’œufs qu’il faut pour constituer une proie 
normale pour nourrir un carnivore. 

Dans le décompte des carnivores nourris, il faut tenir compte de cet 
apport de nourriture supplémentaire. La nouvelle valeur AN' est donnée 
par l'équation : 
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AN'(K) = AN(K) + («? A xzan/POAN) . (0 ANzas /POAS) + (0 ANZ bn 
/POBN) + (0 anzbs /POBS) + (0 A nz C /POC) + (0 axze /POE) 

+ (0anzf/POF) + (WPOG) + ( (ÔaNZTRA\/POAN) 

+ (®AX ZTRAs/POAS) + AN ZTRBn/POBN ) + ( 6AS ZTRHs/POBS ) 

+ (#AN ZTRc/POC + [9\n ZTKe/POE) 

4" (0AX ZTRf/POF) + ( 6\s ZTRG /POG)) . PP)/BE 

— Stades larvaires 

Dans le modèle proposé, l’évolution des stades larvaires est considérée 
comme étant progressive et le terme « jeunes » s’applique à tous ces 
stades. 

A un moment défini (par exemple, K3), les œufs éclosent. En pre¬ 
mière approximation, on va supposer que tous les œufs éclosent en même 
temps. 

Considérons les jeunes carnivores du groupe AN : 

JAN(K3) = ZTRAN(K3) + ZTRAS(K3) + ZAN(K3) + ZAS(K3) 

Pendant la durée de leur développement (jusqu’au temps K4, par 
exemple), les jeunes sont des proies pour les adultes. 

On désigne par GGAMJA(K), un coefficient variable indiquant com¬ 
bien il faut de jeunes pour constituer une proie de dimension normale. 

D’autres variables interviennent également : 

— le diamètre de perception ou d’influence DJA<DA; 

— le coefficient de correction prédateur-proie EEPSJA (les jeunes ne 
peuvent manger toutes les proies qu’ils rencontrent, vu les tailles 
souvent trop importantes des proies); 

— la vitesse de déplacement VJA<VA. 

A partir du moment K4, les jeunes sont adultes. A ce moment, on ne 
les considère plus comme proies et les valeurs de DJA, VJA et EEPSJA 
sont identiques à celles des adultes. Ce passage du stade jeune au 
stade adulte est progressif : DJA, VJA, EEPSJA et GGAMJA varient, 
en fonction du temps, de façon, en première approximation, linéaire. 

La durée du phénomène est K4 — K3. On a les équations suivantes : 

(DA - DJA(K3) ) 

DJA(K) = DA +-. K 

(K4 - K3) 



47, 19 


DANS 2 ZONES CLIMATIQUEMENT DIFFÉRENTES 


31 


(GGAMJA(K3) - 1) (K4 - K3 - K) 

GGAMJA(K) =-+ 1 

K4 - K3 

(EEPSJA(K) - EEPSJA(K3) ) 

EEPSJA(K) = EEPSJA(K3) +-. K 

K4 - K3 

L'équation donnant la valeur de VJA en fonction du temps est de 
même forme que celle de DJA. 

Comme précédemment défini pour les animaux adultes, les jeunes ont 
le choix entre 2 milieux distincts, de caractéristiques écologiques diffé¬ 
rentes. Leur distribution dans l’un ou l'autre milieu est également sous 
la dépendance de leur préférendums. Un système d’équations analogue 
à celui défini dans le cas des adultes peut donc être utilisé, avec des pré¬ 
férendums identiques ou différents. 

De même, le nombre d'individus vivant dans le milieu est sous la 
dépendance des mêmes facteurs, et notamment, de la prédation. Les 
équations utilisées dans ces calculs sont du même type que celles utilisées 
pour les adultes avec quelques additions supplémentaires. En effet, les 
jeunes carnivores sont des proies pour les adultes et leurs caractéristiques 
varient en fonction de leur développement. 

A un certain moment (K4, par exemple), les jeunes sont adultes et 
sont incorporés aux populations adultes existantes. 

Le cycle précédemment défini peut alors recommencer. 


II. PROGRAMMATION DES EQUATIONS 
1. Méthodes 

La programmation des équations a été effectuée en langage machine 
FORTRAN IV. La simulation du modèle écologique est constituée de 
35 sous-programmes, détaillés plus loin. Cette subdivision est avantageuse 
lors de la compilation et lors de la mise au point. En outre, elle permet 
d'apporter des modifications ou des améliorations aux sous-programmes 
sans altérer l'ensemble. Les relations entre les paramètres ne sont, en 
effet, pas considérées comme définitives et seront modifiées suivant les 
besoins. 

Les données, les résultats intermédiaires et définitifs sont conservés 
à travers tout le programme dans des zones de mémoires communes 
(« COMMON »). Les valeurs de chaque population dans la période 
étudiée est. de plus, conservée sur bande pour utilisation ultérieure. C'est 
cette bande qui est employée pour la représentation graphique des résul¬ 
tats en utilisant la table traçante BENSON stimulée par les sous-pro- 
grammes appropriés. 
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2, Description des programmes utilisés. 

Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 

1 PP1 710 1. Délimite les noms et les vecteurs com¬ 

muns ainsi que les variables communes. 

2. Coordonne les n travaux de la phase 
PHI en appelant les sous-routines 
(SUBR.) : 

- EFFAC 

- COSSIN (1 seul appel) 

— MOU (autant d'appels que de po¬ 
pulation considérées) 

3. Appelle la phase PH2. 

1 EFFA1 3F5 1. Remet les valeurs de toutes les mémoi¬ 

res à zéro. 

1 COSSINl 3F4 1. Calcule 100 valeurs particulières d'un 

sinus variant de 0 à ir par pas de 

r/100 et les range dans SINN. 

2. Calcule 100 valeurs particulières d'un 
cosinus variant de — r/2 à r/2 par pas 
de r/100 et les range dans COSS. 

3. Imprime des messages signalant que 
ces travaux sont exécutés. 

1 MOU1 5AO 1. Effectue la lecture des données concer¬ 

nant les préférendums des animaux 
étudiés. 

2. Appelle successivement des sous-rou¬ 
tines : 

- CALCL 

- GAMT 
_ GAMH 

1 CALCL1 480 1. Imprime des messages signalant que 

les préférendums vont être calculés. 

2. Appelle la sous-routine PREF 2 fois. 


Travaux exécutés 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


1 PREFl 


1 GAMTl 


1 TTREL1 


1 GAMHl 


1 RRHELl 


1 RRPA1 


478 


3AO 


328 


360 


36C 


390 


1. Calcule les coefficients de répartition 
en fct de la température pour une espè¬ 
ce donnée en appellant RRPA. 

2. Calcule 100 valeurs particulières du 
coefficient de répartition-température 
et les imprime par pas de 10. 

3. Calcule 100 valeurs particulières du 
coefficient de répartition-humidité et 
les imprime par pas de 10. 

1. Calcule les valeurs des gammas pour la 
température en appelant la subroutine 
TTREL. 

2. Calcule 100 valeurs particulières et les 
imprime par pas de 10. 

1. Imprime un message indiquant le pas¬ 
sage dans ce sous-programme. 

2. Calcule le pas des abscisses. 

3. Calcule 100 valeurs d'abscisse. 

4. Calcule 100 valeurs de gammas-tem- 
pérature. 

1. Calcule 100 valeurs de gammas pour 
l’humidité en appelant la sous-routine 
RRHEL. 

2. Calcule 100 valeurs particulières et les 
imprime par pas de 10. 

1. Imprime un message indiquant le pas¬ 
sage dans ce sous-programme. 

2. Calcule 100 valeurs du gamma-humi¬ 
dité. 

1. Imprime un message indiquant le pas¬ 
sage dans ce sous-programme. 

2. Calcule le pas des abscisses. 

3. Calcule 100 valeurs particulières des 
abscisses. 
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Nb bytes 

Phase Nom Subr. (syst. „ 

N o hexa _ Travaux executes 

décimal) 

4. Calcule la valeur absolue des coeffi¬ 
cients de répartition. 

5. Calcule les coefficients en pourcentage. 

2 PP2 1CCO 1. Effectue la lecture des données sui¬ 

vantes : 

— indices temps 

— coefficients divers (surface, popu¬ 
lation, etc.) 

— valeurs des populations intiales 

— paramètres biologiques 

2. Imprime les valeurs des données. 

3. Calcule la variation de vitesse en fct 
d'un coefficient de sous-alimentation. 

4. Calcule les coefficients de mortalité et 
les imprime. 

5. Calcule 10.000 valeurs intervenant 
dans une intégrale de la sous-routine 
TETA. 

6. Imprime 10 valeurs de S intervenant 
dans TETA. 

7. Imprime les données. 

8. Appelle la sous-routine TEMP. 

Les activités 9, 10. 11 suivantes sont appe¬ 
lées autant de fois que demandées par la 
valeur de K5. 

9. Appelle la sous-routine DISTRIB. 

10. Appelle la sous-routine VARIAT avec 
les arguments intervenant dans la zone 
extérieure. 

11. Appelle la sous-routine VARIAT avec 
les arguments intervenant dans la zone 
profonde. 

12. Appelle la phase PH3. 

2 TEMP1 2DC 1. Calcule 30 valeurs particulières de la 

température et de l’humidité dans cha¬ 
cune des 2 zones. 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


2 DSTRIBl A8C 

2 DISTRI 7B8 


2 RPARTI1 AA8 


2 GLOBAL1 968 


2. Imprime ces 120 valeurs. 

NB : ce programme est à modifier lorsque 
l’on veut modifier les valeurs des condi¬ 
tions physiques. 

1. Appelle la sous-routine DISTR au¬ 
tant de fois qu’il y a de populations 
expérimentales. 

1. Vérifie si les conditions physiques uti¬ 
lisées ne sont pas léthales. Si oui, un 
message est envoyé et le programme 
est arrêté. 

2. Appelle la sous-routine RPARTI avec 
les arguments faisant intervenir la tem¬ 
pérature. 

3. Appelle la sous-routine RPARTI avec 
les arguments faisant intervenir l’hu¬ 
midité. 

4. Appelle la sous-routine GLOBAL. 

5. Effectue les calculs donnant les va¬ 
leurs de distribution de la population 
considérée, dans les 2 zones. 

6. Imprime les valeurs trouvées en 5). 

1. Calcule les coefficients de répartition 
RTE et RTP en fonction des condi¬ 
tions physiques données. 

2. Imprime les valeurs de ces coefficients. 

1. Imprime la valeur des paramètres phy¬ 
siques au moment considéré. 

2. Vérifie si la valeur trouvée n'est pas 
léthale. Si oui, le travail est arrêté. 

3. Calcule la valeur des 2 gammas dans 
ces conditions pour chacune des 2 
zones. 

4. Calcule les valeurs globalisées de ces 
4 gammas et les imprime. 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


2 VARIATl C70 


2 VAR1A21 7FC 


2 TETAI 300 

2 RESULTl 2FC 

2 SATUR1 554 

3 PP3 232C 


1. Appelle la sous-routine VARIA2 au¬ 
tant de fois qu'il y a de populations 
de carnivores. 

2, Calcule la valeur des différentes popu¬ 
lations considérées. 

1. Calcule les variations de populations 
en fct du nombre de rencontres, en uti¬ 
lisant la fonction TETA, 

2. Vérifie si les carnivores mangent plus 
d'une proie par unité de temps. Si oui, 
voir 3 ). 

3. Calcule le nombre de proies mangées 
par unité de temps. 

4. Calcule en appelant la sous-routine 
RESULT. 

5. Vérifie si le nb de proies mangées ne 
dépasse pas une certaine valeur maxi¬ 
mum PRGX, Si oui, voir 6. 

6. Appelle la sous-routine SATURA. 

7. Vérifie si la valeur de la population 
étudiée n'est pas nulle ou négative. 

1. Calcule la valeur du nb de rencontre 
entre 2 populations en considérant di¬ 
vers facteurs biologiques : vitesse, dia¬ 
mètre de perception, densité. 

1. Imprime les valeurs de PRMGAN au 
moment considéré. 

1. Imprime un message indiquant que 
l'appétit des carnivores est saturé. 

2. Effectue les corrections nécessaires. 

3. Imprime la valeur de PRGXAN. 

1. Calcule les valeurs des vitesses mini¬ 
males des jeunes. 

2. Calcule les valeurs des diamètres mini- 
ma des jeunes. 
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Nb bytes 

N* (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


3 CFFl 
3 ACCOU1 


3. En utilisant la sous-routine fonction 
CFF, calcule les coefficients : 

— nb. de jeunes pour constituer une 
proie 

— variation de vitesse des jeunes 

— coefficient de correction prédateur 
proie pour les jeunes carnivores. 

— variation du diamètre des jeunes 
Ces calculs sont effectués pendant la pé¬ 
riode juvénile fixée. 

A. Imprime les valeurs des populations 
pendant toute la période expérimentale. 

5. Appelle la sous-routine ACCOUP 2 
fois en utilisant les arguments de cha¬ 
que zone. 

6. Appelle la sous-routine OEUFS 
comme en 5. 

7. Imprime les nouvelles valeurs de po¬ 
pulations. 

8. Imprime les valeurs des populations 
d’œufs chez toutes les espèces consi¬ 
dérées. pour chacune des 2 zones. 

9. Calculent le nb total de jeunes qui 
éclosent. 

10. Imprime ces valeurs pour chacune des 
2 zones. 


220 1. Imprime les valeurs d'un coefficient 

s'appliquant aux jeunes. 

1EB8 1. Vérifie la longueur des vecteurs utili¬ 
sés. En cas de dépassement de capa¬ 
cité, le programme est arrêté et un 
message est envoyé. 

2. Imprime un message signalant le pas¬ 
sage dans ce sous-programme. 

3. Calcule le nb de mâles et de femelles 
dans chacune des populations. 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


3 OEUF1 


4 PP4 


3D74 


D98 


3. Calcule le nb de mâles et de femelles 
dans chacune des populations. 

4. Imprime ces valeurs. 

5. Calcule le nb d’accouplements pen¬ 
dant la période considérée en utilisant 
le sous-programme TETA. 

6. Calcule le nb de femelles fécondées. 

7. Imprime le nb d’accouplements pour 
chaque population. 

8. Calcule la proportion de femelles fé¬ 
condées. 

9. Calcule le nb d’œufs pondus. 

1. Imprime un message indiquant que ce 
sous-programme est utilisé. 

2. Vérifie la longueur des vecteurs utili¬ 
sés. En cas de dépassement de capa¬ 
cité, le programme est arrêté et un 
message est envoyé. 

3. Calcule le nb d’œufs en fct du temps 
en considérant le nb d’attaques de 
prédateurs. Utilise le sous-programme 
TETA, vérifie si le nb d'œufs ne 
devient pas négatif. 

4. Imprime les valeurs des indices temps. 

5. Calcule le nb de carnivores nourris 
par l’absorption d’œufs. 

6. Imprime les résultats partiels. 

1. Calcule la distribution des jeunes dans 
les 2 zones en appelant la sous-routine 
DSTRIJ. 

2. Calcule les variations du nb de jeunes 
en fct du temps en appelant 2 fois 
JARIAT affecté des arguments de 
chacune des 2 zones. 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


4 DSTRIJ1 


4 JARIAl 


5 PP5 


3. Calcule le nb total de jeunes. 

N.B. : 1, 2, 3 sont calculés autant de fois 
que le nécessite les valeurs de l’indice 
temps K4, K5. 

4. Imprime les valeurs de données à titre 
de vérification. 

5. Appelle PH5. 

A8C 1. Travaille de façon analogue à DISTR 
en appelant les mêmes sous-program¬ 
mes : 

DISTR 

RPARTI 

GLOBAL 


2FDC 1. Vérifie la longueur des vecteurs utili¬ 
sés et imprime un message en cas de 
dépassement de capacité avec arrêt 
éventuel du travail. 

2. Calcule les variations de chaque popu¬ 
lation de jeunes en utilisant le sous- 
programme TETA. En cas d’appari¬ 
tion de valeurs nulles ou négatives, le 
calcul est arrêté et le contrôle est ren¬ 
voyé au programme principal PP3. 

3. La valeur des populations des jeunes 
au temps K5 sont ajoutées aux popu¬ 
lation d'adultes. 


1CCO 1. Exécute les mêmes calculs que PP2 
mais dans l’intervalle de temps K5. Il, 
en utilisant les mêmes sous-routines, à 
l'exclusion de TEMP : 

- DISTRB 

- DISTR 

- RPARTI 

- GLOBAL 

- VARIAT 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


6 PP6 


6 BANDE1 


6 AMAX1 

7 PP7 


7 BENSON 


- VARIA2 

- TETA 

- RESULT 

- SATURA 
2. Appelle PH6. 


334 1. Remet la bande réservée sur l'unité 

n" 2 au départ. 

2. Appelle 2 fois la sous-routine BANDE 
pour chacune des 2 zones, avec les 
arguments correspondants. 

3. Appelle PH7, 


16DO 


1. Imprime un message indiquant que ce 
sous-programme est utilisé. 

2. Remet certaines valeurs de populations 
à 0 pour simplifier le traçage. 

3. Ecrit les valeurs des populations sur 
la bande réservée n° 2. 

4. Calcule la valeur maximum de n'im¬ 
porte quelle population expérimentale 
en utilisant AMAX1. 

5. Imprime cette valeur sur la bande ré¬ 
servée n" 2 et l'imprime. 


Sous-routine de biliothèque. 

318 1. Remet la bande réservée n" 2 au dé¬ 

part. 

2. Calcule la valeur maximum de l’or¬ 
donnée des 2 graphiques AMAX. 

3. Appelle la sous-routine BENSON 
pour chacune des 2 zones. 

4. Termine le travail par CALL EXIT. 

1A64 1. Calcule certaines valeurs pour le cal¬ 

cul des coordonnées. 
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Nb bytes 

N° (syst. 

Phase Nom Subr. hexa¬ 
décimal) 


Travaux exécutés 


2. Lit les valeurs des populations étudiées 
à partir de la bande se trouvant sur 
l'unité n° 2. 

3. Calcule les valeurs des abscisses du 
graphique. 

4. Appelle la sous-routine NUTRAC 
avec le n° du tracé. 

5. Appelle les sous-routines de biblio¬ 
thèque AXES, BENINT et PRINT. 

6. Appelle TRACE. 

7. Imprime un message indiquant que le 
travail BENSON est terminé. 

7 TRACEl 1 324 1. Calcule les 2 coefficients intervenant 

dans les coordonnées. 

2 Appelle successivement les sous-rou¬ 
tines : 

- COURBE 

- BENINT 

- PRINT. 

Remarques complémentaires 

1. La place en nbX de bytes occupées par les programmes principaux ne 
comprend pas la place utilisée par les COMMON (valeurs commu¬ 
nes) et qui est de Fl FC. 

2. Les noms des sous-routines repris dans la colonne « NOM SUB. » 
sont les noms utilisés dans les cartes contrôle et désignent donc les 
noms des MODULES binaires. Dans la rubrique « Travaux exécu¬ 
tés », les noms des sous-routines désignent les noms utilisés dans les 
appels CALL correspondants. 



42 


R. R. TERCAFS. — SIMULATION DUN SYSTÈME ÉCOLOGIQUE 


47, 19 


3. Assemblage des programmes 

Les programmes peuvent être assemblés de façon différente suivant 
le système d’édition disponible. Nous avons ainsi successivement utilisé 
la procédure PS 44 sur IBM 360/44 (séparation en phases) et la procé¬ 
dure OS sur IBM 360/65 (additions successives avec programme principal 
appelant; « root phase overlay » ). 


a. Edition en PS 44 
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b) Edition en OS 


EFFAC 

COSSIN 



PP7 

LECT 

BENSON 


4. Programme FORTRAN 

Les documents relatifs au programme FORTRAN, faisant l'objet d’une 
programmation sur ordinateur, peuvent être obtenus sur demande au 
Patrimoine de l'Institut royal des Sciences naturelles de Belgique. 


5. Abréviations utilisées 


Symbole 

Longueur 

vecteur 


Signification 

S 

101 j 

Valeur intervenant dans le calcul du nombre de 

A 

102 j 

rencontres. 


SINN 

100 

Contient 100 

valeurs particulières d'un sinus 


variant de 0 à ?r par pas de tt/100. 
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Symbole 


Longueur 

vecteur 


Signification 


coss 

100 

Concient 100 valeurs particulières d'un cosinus 
variant de 0 à tt par pas de 7r/100. 

ZERP 

30 

Coefficient de répartition de la : 

zone extérieure. 

ZPRP 

30 

Coefficient de répartition de la 

zone profonde. 

ANVLO 

1 

Valeur minimum où le gamma 
animaux AN existe. 

thermique 

des 

ANVLT 

1 

Valeur maximum où le gamma 
animaux AN existe. 

thermique 

des 

BNVLO 

1 

Valeur minimum où le gamma 
animaux BN existe. 

thermique 

des 

BNVLT 

1 

Valeur maximum où le gamma 
animaux BN existe. 

thermique 

des 

ASVLO 

1 

Valeur minimum où le gamma 
animaux AS existe. 

thermique des 

ASVLT 

1 

Valeur maximum où le gamma 
animaux AS existe. 

thermique 

des 

BSVLO 

1 

Valeur minimum où le gamma 
animaux BS existe. 

thermique 

des 

BSVLT 

1 

Valeur maximum où le gamma 
animaux BS existe. 

thermique 

des 

CCVLO 

1 

Valeur minimum où le gamma 
animaux C existe. 

thermique 

des 

CCVLT 

1 

Valeur maximum où le gamma 
animaux C existe. 

thermique 

des 

EEVLO 

1 

Valeur minimum où le gamma 
animaux E existe. 

thermique 

des 

EEVLT 

1 

Valeur maximum où le gamma 
animaux E existe. 

thermique des 

FFVLO 

1 

Valeur minimum où le gamma 
animaux F existe. 

thermique 

des 

FFVLT 

1 

Valeur maximum où le gamma 
animaux F existe. 

thermique 

des 

GGVLO 

1 

Valeur minimum où le gamma 

thermique 

des 


animaux G existe. 
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Symbole 

GGVLT 

ANLOT 

ANLTT 

ANVLMT 

ANLOH 

ANLTH 

ANVLMH 

ANRLMI 

ANRLMX 

ANRHO 

ANRELM 

ANREMI 

ASLOT 

ASLTT 

ASVLMT 

ASLOH 


Longueur 

vecteur 


Signification 


1 Valeur maximum où le gamma thermique des 
animaux G existe. 

1 Valeur minimum où le thermopréférendum des 
animaux AN existe. 

1 Valeur maximum où le thermopréférendum des 
animaux AN existe. 

1 Valeur maximum de l'intensité du thermopréfé¬ 
rendum des animaux AN. 

1 Valeur minimum où l'hygropréférendum des 

animaux AN existe. 

1 Valeur maximum où l’hygropréférendum des 

animaux AN existe. 

1 Valeur de l'intensité maximum de l'hygropréfé- 
rendum des animaux AN. 

1 Valeur de l’intensité minimum du gamma ther¬ 
mique des animaux AN. 

1 Valeur de l’intensité maximum du gamma ther¬ 
mique des animaux AN. 

1 Valeur du point d’inflexion du gamma hygro¬ 
métrique des animaux AN. 

1 Valeur de l'intensité maximum du gamma hygro¬ 
métrique des animaux AN. 

1 Valeur de l'intensité minimum du gamma hygro¬ 
métrique des animaux AN. 

1 Valeur minimum où le thermopréférendum des 
animaux AS existe. 

1 Valeur maximum où le thermopréférendum des 
animaux AS existe. 

1 Valeur de l'intensité maximum du thermopréfé¬ 
rendum des animaux AS. 

1 Valeur minimum où l'hygropréférendum des ani¬ 
maux AS existe. 
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Symbole 

ASLTH 

ASVLMH 

ASRLMI 

ASRLMX 

ASRHO 

ASRELM 

ASREMI 

B N LOT 

BNLTT 

BNVLMT 

BNLOH 

BNLTH 

BNVLMH 

BNRLMI 

BNRLMX 

BNRHO 

BNRELM 


Longueur 

vecteur 


Signification 


1 Valeur maximum où l’hygropréférendum des 
animaux AS existe. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux AS. 

1 Valeur de l'intensité minimum du gamma ther¬ 
mique des animaux AS. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux AS. 

1 Valeur du point d'inflexion du gamma hygro¬ 
métrique des animaux AS. 

1 Valeur de l’intensité maximum du gamma hygro¬ 
métrique des animaux AS. 

1 Valeur de l'intensité minimum du gamma hygro¬ 
métrique des animaux AS. 

1 Valeur minimum où le thermopréférendum des 
animaux BN existe. 

1 Valeur maximum où le thermopréférendum des 
animaux BN existe. 

1 Valeur maximum de l’intensité du thermopréfé¬ 
rendum des animaux BN. 

1 Valeur minimum où l'hygropréférendum des 
animaux BN existe. 

1 Valeur maximum où l'hygropréférendum des 
animaux BN existe. 

1 Valeur de l'intensité maximum de l'hygropréfé- 
rendum des animaux BN. 

1 Valeur de l'intensité minimum du gamma ther¬ 
mique pour les animaux BN. 

1 Valeur de l’intensité maximum du gamma ther¬ 
mique des animaux BN. 

1 Valeur du point d'inflexion du gamma hygro¬ 
métrique des animaux BN. 

1 Valeur de l’intensité maximum du gamma hygro¬ 
métrique des animaux BN. 
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Symbole 

BNREMI 

BSLOT 

BSLTT 

BSVLMT 

BSLOH 

BSLTH 

BSVLMH 

BS R LM I 

BSRLMX 

BSRHO 

BSRELM 

BSREMI 

CCLOT 

CCLTT 

CCVLMT 

CCLOLI 

CCLTH 


Longueur 

vecteur 


Signification 


1 Valeur de l’intensité minimum du gamma hygro¬ 
métrique des animaux BN. 

1 Valeur minimum où le thermopréférendum des 
animaux BS existe. 

1 Valeur maximum où le thermopréférendum des 
animaux BS existe. 

1 Valeur maximum du thermopréférendum des 
animaux BS. 

1 Valeur minimum où l’hygropréférendum des ani¬ 
maux BS existe. 

1 Valeur maximum où l'hygropréférendum des 
animaux BS existe. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux BS. 

1 Valeur de l'intensité minimum du gamma ther¬ 
mique des animaux BS. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux BS. 

1 Valeur du point d’inflexion du gamma hygro¬ 
métrique des animaux BS. 

1 Valeur de l'intensité maximum du gamma hygro¬ 
métrique des animaux BS. 

1 Valeur de l'intensité minimum du gamma hygro¬ 
métrique des animaux BS. 

1 Valeur minimum où le thermopréférendum des 
animaux C existe. 

1 Valeur maximum où le thermopréférendum des 
animaux C existe. 

1 Valeur maximum de l'intensité du thermopréfé¬ 
rendum des animaux C. 

1 Valeur minimum où l’hygropréférendum des 
animaux C existe. 

1 Valeur maximum où l’hygropréférendum des 
animaux C existe. 
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Symbole 

CCVLMH 

CCRLMI 

CCRLMX 

CCRHO 

CCRELM 

CCREMI 

EELOT 

EELTT 

EEVLMT 

EELOH 

EELTH 

EEVLMH 

EERLMI 

EERLMX 

EERHO 

EERELM 

EEREMI 


Longueur 

vecteur 


Signification 


1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux C. 

1 Valeur de l'intensité minimum du gamma ther¬ 
mique des animaux C. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux C. 

1 Valeur du point d'inflexion du gamma hygro¬ 
métrique des animaux C. 

1 Valeur de l'intensité maximum du gamma hygro¬ 
métrique des animaux C. 

1 Valeur de l'intensité minimum du gamma hygro¬ 
métrique des animaux C. 

1 Valeur minimum où le thermopréférendum des 
animaux E existe. 

1 Valeur maximum où le thermopréférendum des 
animaux E existe. 

1 Valeur maximum de l'intensité du thermopréfé¬ 
rendum des animaux E. 

1 Valeur minimum où l'hygropréférendum des 
animaux E existe. 

1 Valeur maximum où l’hygropréférendum des 
animaux E existe. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux E. 

1 Valeur de l’intensité minimum du gamma ther¬ 
mique des animaux E. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux E. 

1 Valeur du point d'inflexion du gamma hygro¬ 
métrique des animaux E. 

1 Valeur de l'intensité maximum du gamma hygro¬ 
métrique des animaux E. 

1 Valeur de l'intensité minimum du gamma hygro¬ 
métrique des animaux E. 
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Symbole 

F F LOT 

FFLTT 

FFVLMT 

FFLOH 

FFLTH 

FFVLMH 

FFRLMI 

FFRLMX 

FFRHO 

FFRELM 

FFREMI 

GGLOT 

GGLTT 

GGVLMT 

GGLOH 

GGLTH 

GGVLMH 


Longueur 

vecteur 


Signification 


1 Valeur minimum où le thermopréférendum des 
animaux F existe. 

1 Valeur maximum où le thermopréférendum des 
animaux F existe. 

1 Valeur maximum de l’intensité du thermopréfé¬ 
rendum des animaux F. 

1 Valeur minimum où l'hygropréférendum des 
animaux E existe. 

1 Valeur maximum où l'hygropréférendum des 
animaux E existe. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux F. 

1 Valeur de l'intensité minimum du gamma ther¬ 
mique des animaux F. 

1 Valeur de l’intensité maximum du gamma ther¬ 
mique des animaux F. 

1 Valeur du point d’inflexion du gamma hygro¬ 
métrique des animaux F. 

1 Valeur de l’intensité maximum du gamma hygro¬ 
métrique des animaux F. 

1 Valeur de l’intensité minimum du gamma hygro¬ 
métrique des animaux F. 

1 Valeur minimum où le thermopréférendum des 
animaux G existe. 

1 Valeur maximum où le thermopréférendum des 
animaux G existe. 

1 Valeur maximum de l’intensité du thermopréfé¬ 
rendum des animaux G. 

1 Valeur minimum où l'hygropréférendum des 
animaux G existe. 

1 Valeur maximum où l’hygropréférendum des 
animaux G existe. 

1 Valeur de l'intensité maximum du gamma ther¬ 
mique des animaux G. 
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Symbole 


Longueur 

vecteur 


Signification 


GGRLMI 

1 Valeur de l'intensité minimum du gamma ther¬ 

mique des animaux G. 

GGRLMX 

1 Valeur de l'intensité maximum du gamma ther¬ 

mique des animaux G. 

GGRHO 

1 Valeur du point d’inflexion du gamma hygro¬ 

métrique des animaux G. 

GGRELM 

1 Valeur de l'intensité maximum du gamma hygro¬ 

métrique des animaux G. 

GGREMI 

1 Valeur de l’intensité minimum du gamma 

hygrométrique des animaux G. 

SALA 

1 Coefficient de sous-alimentation affectant la 

vitesse des animaux AN. 

SALB 

1 Coefficient de sous-alimentation affectant la 

vitesse des animaux BN. 

DA 

1 Diamètre de perception ou d’influence des ani¬ 

maux du groupe A. 

DB 

1 Diamètre de perception ou d'influence des ani¬ 

maux du groupe B. 

DC 

1 Diamètre de perception ou d'influence des ani¬ 

maux C. 

DE 

1 Diamètre de perception ou d'influence des ani¬ 

maux E. 

DF 

1 Diamètre de perception ou d’influence des ani¬ 

maux du groupe F. 

DG 

1 Diamètre de perception ou d'influence des ani¬ 

maux du groupe G. 

VA 

1 Vitesse de déplacement des animaux du 

groupe A. 

VB 

1 Vitesse de déplacement des animaux du 

groupe B. 

VC 

1 Vitesse de déplacement des animaux C. 

VE 

1 Vitesse de déplacement des animaux E. 

VF 

1 Vitesse de déplacement des animaux F. 
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Symbole 


Longueur 

vecteur 


Signification 


VG 1 Vitesse de déplacement des animaux G. 

NA 1 « Sex-ratio » chez les animaux du groupe A. 

NB 1 « Sex-ratio » chez les animaux du groupe B. 

NC 1 « Sex-ratio » des animaux du groupe C. 

NE 1 « Sex-ratio » chez les animaux du groupe E. 

N F 1 « Sex-ratio » chez les animaux du groupe F. 

NG 1 « Sex-ratio » chez les animaux du groupe G. 

TRA 1 Proportion d’espèces creusant des terriers chez 

les animaux du groupe A. 

TRB 1 Proportion d'espèces creusant des terriers chez 

les animaux du groupe B. 

TRC 1 Proportion d'espèces creusant des terriers chez 

les animaux du groupe C. 

TRE 1 Proportion d'espèces creusant des terriers chez 

les animaux E. 

TRF 1 Proportion d’espèces creusant des terriers chez 

les animaux F. 

TRG 1 Proportion d'espèces creusant des terriers chez 

les animaux G. 

OTRAN 1 Nombres d’œufs pondus par les espèces à ter¬ 

rier chez les animaux AN. 

OTRAS 1 Nombre d'œufs pondus par les espèces à ter¬ 

rier chez les animaux AS. 

OTRBN 1 Nombre d'œufs pondus par les espèces à ter¬ 

rier chez les animaux BN, 

OTRBS 1 Nombre d'œufs pondus par les espèces à ter¬ 

rier chez les animaux BS. 

OTRC 1 Nombre d'œufs pondus par les espèces à ter¬ 

rier chez les animaux C. 

OTRE 1 Nombre d'œufs pondus par les espèces à ter¬ 

rier chez les animaux E. 
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Symbole 


Longueur 

vecteur 


Signification 


OTRF 

OTRG 

OAN 

OAS 

OBS 

OBN 

OC 

OE 

OF 

OG 

DZA 

DZB 

DZC 

DZE 

DZF 

DZG 

POAN 

POAS 

POBS 

POBN 

POC 

POE 

POF 

POG 

BE 


1 Nombre d'œufs pondus par les espèces à ter¬ 
rier chez les animaux F. 

1 Nombre d'œufs pondus par les espèces à ter¬ 
rier chez les animaux G. 

1 Nombres d'œufs pondus par les animaux AN. 

1 Nombre d’œufs pondus par les animaux AS, 

1 Nombre d’œufs pondus par les animaux BS. 

1 Nombre d'œufs pondus par les animaux BN, 

1 Nombre d'œufs pondus par les animaux C. 

1 Nombre d'œufs pondus par les animaux E. 

1 Nombre d'œufs pondus par les animaux F. 

1 Nombre d'œufs pondus par les animaux G. 

1 Diamètre des couvées des animaux du groupe A, 

1 Diamètre des couvées des animaux du groupe B. 

1 Diamètre des couvées des animaux C. 

1 Diamètre des couvées des animaux E. 

1 Diamètre des couvées des animaux F. 

1 Diamètre des couvées des animaux G, 

1 Nombre d'œufs dans une ponte d'animaux AN. 

1 Nombre d’œufs dans une ponte d'animaux AS. 

1 Nombre d'œufs dans une ponte d'animaux BS. 

1 Nombre d'œufs dans une ponte d'animaux BN. 

1 Nombre d’œufs dans une ponte d’animaux C. 

1 Nombre d'œufs dans une ponte d’animaux E. 

1 Nombre d’œufs dans une ponte d'animaux F. 

1 Nombre d'œufs dans une ponte d'animaux G. 

1 Nombre d'œufs qu'il faut pour constituer une 

proie de dimension normale. 
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Symbole 


Longueur 

vecteur 


Signification 


PP 

GAMJA 

GAMJB 

GAMJC 

GAMJE 

GAMJF 

GAMJC 

IJ 

EPSJA 

EPSJB 

II 

Kl 

K2 

K3 

K4 

K5 

Y 

X 

U1 

PRGXAN 


1 Coefficient de protection dû au terrier. 

1 Nombre de jeunes de l'espèce A nécessaire pour 

constituer une proie normale. 

1 Nombre de jeunes du groupe B nécessaire pour 
constituer une proie normale. 

1 Nombre de jeunes du groupe C nécessaire pour 
constituer une proie normale. 

1 Nombre de jeunes du groupe E nécessaire pour 
constituer une proie normale. 

1 Nombre de jeunes du groupe F nécessaire pour 
constituer une proie normale. 

1 Nombre de jeunes du groupe G nécessaire pour 
constituer une proie normale. 

1 Coefficient affectant la vitesse des jeunes. 

1 Coefficient de protection prédateur-proie pour 
les jeunes du groupe A. 

1 Coefficient de protection prédateur-proie pour 
les jeunes du groupe B. 

1 Valeur maximum de l’indice temps. 

1 Début de la période d'accouplement. 

1 Fin de la période d’accouplement, début de la 
gestation. 

1 Fin de la gestation, ponte des œufs. 

1 Moment de l'éclosion, apparition des jeunes. 

1 Fin de la période juvénile. 

1 Coefficient de mortalité. 

1 Coefficient de densité de peuplement. 

1 Surface de la zone étudiée. 

1 Nombre de proies maximum mangeable par un 

carnivore AN pendant une unité de temps. 
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Symbole 

Longueur 

vecteur 

Signification 

PRGXAS 

1 

Nombre de proies maximum mangeable par un 
carnivore AS pendant une unité de temps. 

PRGXBN 

1 

Nombre de proies maximum mangeable par un 
carnivore BN pendant une unité de temps. 

PRGXBS 

1 

Nombre de proies maximum mangeable par un 
carnivore BS pendant une unité de temps. 

EZTRAN 

30 

Nombre d'œufs des espèces à terrier du groupe 
AN dans la zone extérieure. 

EZTRBN 

30 

Nombre d'œufs des espèces à terrier du groupe 
BN dans la zone extérieure. 

EZTRAS 

30 

Nombre d’œufs des espèces à terrier du groupe 
AS dans la zone extérieure. 

EZTRBS 

30 

Nombre d'œufs des espèces à terrier du groupe 
BS dans la zone extérieure. 

EZTRC 

30 

Nombre d'œufs des espèces à terrier du groupe 
C dans la zone extérieure. 

EZTRE 

30 

Nombre d'œufs des espèces à terrier du groupe 
E dans la zone extérieure. 

EZTRF 

30 

Nombre d’œufs des espèces à terrier du groupe 
F dans la zone extérieure. 

EZTRG 

30 

Nombre d'œufs des espèces à terrier du groupe 
G dans la zone extérieure. 

PZTRAN 

30 

Nombre d’œufs des espèces à terrier du groupe 
AN dans la zone profonde. 

PZTRAS 

30 

Nombre d'œufs des espèces à terrier du groupe 
AS dans la zone profonde. 

PZTRBN 

30 

Nombre d’œufs des espèces à terrier du groupe 
BN dans la zone profonde. 

PZTRBS 

30 

Nombre d’œufs des espèces à terrier du groupe 
BS dans la zone profonde. 

PZTRC 

30 

Nombre d’œufs des espèces à terrier du groupe 
C dans la zone profonde. 
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Symbole 

Longueur 

vecteur 

Signification 

PZTRE 

30 

Nombre d'œufs des espèces à terrier du groupe 
E dans la zone profonde. 

PZTRF 

30 

Nombre d'œufs des espèces à terrier du groupe 
F dans la zone profonde. 

PZTRG 

30 

Nombre d'œufs des espèces à terrier du groupe 
G dans la zone profonde. 

EZAN 

30 

Nombre d'œufs des animaux AN dans la zone 
extérieure. 

EZAS 

30 

Nombre d'œufs des animaux AS dans la zone 
extérieure. 

EZBN 

30 

Nombre d'œufs des animaux BN dans la zone 
extérieure. 

EZBS 

30 

Nombre d’œufs des animaux BS dans la zone 
extérieure. 

EZC 

30 

Nombre d'œufs des animaux C dans la zone 
extérieure. 

EZE 

30 

Nombre d'œufs des animaux E dans la zone 
extérieure. 

EZF 

30 

Nombre d’œufs des animaux F dans la zone 
extérieure. 

EZG 

30 

Nombre d'œufs des animaux G dans la zone 
extérieure. 

PZAN 

30 

Nombre d'œufs des animaux AN dans la zone 
profonde. 

PZAS 

30 

Nombre d’œufs des animaux AS dans la zone 
profonde. 

PZBN 

30 

Nombre d’œufs des animaux BN dans la zone 
profonde. 

PZBS 

30 

Nombre d’œufs des animaux BS dans la zone 
profonde. 

PZC 

30 

Nombre d'œufs des animaux C dans la zone 


profonde. 
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Symbole 

Longueur 

vecteur 

Signification 

PZE 

30 

Nombre d'œufs des animaux E dans la zone 
profonde. 

PZF 

30 

Nombre d'œufs des animaux F dans la zone 
profonde. 

PZG 

30 

Nombre d'œufs des animaux G dans la zone 
profonde. 

ZEAN 

30 

Nombre d'animaux AN dans la zone extérieure. 

ZEAS 

30 

Nombre d'animaux AS dans la zone extérieure. 

ZEBS 

30 

Nombre d'animaux BS dans la zone extérieure. 

ZEBN 

30 

Nombre d'animaux BN dans la zone extérieure. 

ZEC 

30 

Nombre d'animaux C dans la zone extérieure. 

ZEF 

30 

Nombre d’animaux F dans la zone extérieure. 

ZEG 

30 

Nombre d’animaux G dans la zone extérieure. 

ZPAN 

30 

Nombre d'animaux AN dans la zone profonde. 

ZPAS 

30 

Nombre d’animaux AS dans la zone profonde. 

zpbn 

30 

Nombre d'animaux BN dans la zone profonde. 

ZPBS 

30 

Nombre d’animaux BS dans la zone profonde. 

ZPC 

30 

Nombre d’animaux C dans la zone profonde. 

ZPE 

30 

Nombre d’animaux E dans la zone profonde. 

ZPF 

30 

Nombre d’animaux F dans la zone profonde. 

ZPG 

30 

Nombre d'animaux G dans la zone profonde. 

ZEANC 

30 

Nombre de rencontres entre les animaux AN 
et C dans la zone extérieure. 

ZEANE 

30 

Nombre de rencontres entre les animaux AN 
et E dans la zone extérieure. 

ZEANF 

30 

Nombre de rencontres entre les animaux AN 
et F dans la zone extérieure. 

ZEANG 

30 

Nombre de rencontres entre les animaux AN 


et G dans la zone extérieure. 
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Symbole 

Longueur 

vecteur 

Signification 



ZPANG 

30 

Nombre de rencontres entre les 
et C dans la zone profonde. 

animaux 

AN 

ZPANE 

30 

Nombre de rencontres entre les 
et E dans la zone profonde. 

animaux 

AN 

ZPANF 

30 

Nombre de rencontres entre les 
et F dans la zone profonde. 

animaux 

AN 

ZPANG 

30 

Nombre de rencontres entre les 
et G dans la zone profonde. 

animaux 

AN 

ZEASC 

30 

Nombre de rencontres entre les 
et C dans la zone extérieure. 

animaux 

AS 

ZEASE 

30 

Nombre de rencontres entre les 
et E dans la zone extérieure. 

animaux 

AS 

ZEASF 

30 

Nombre de rencontres entre les 
et F dans la zone extérieure. 

animaux 

AS 

ZEASG 

30 

Nombre de rencontres entre les 
et G dans la zone extérieure. 

animaux 

AS 

ZPASC 

30 

Nombre de rencontres entre les 
et C dans la zone profonde. 

animaux 

AS 

ZPASE 

30 

Nombre de rencontres entre les 
et E dans la zone profonde. 

animaux 

AS 

ZPASF 

30 

Nombre de rencontres entre les 
et F dans la zone profonde. 

animaux 

AS 

ZPASG 

30 

Nombre de rencontres entre les 
et G dans la zone profonde. 

animaux 

AS 

ZPBNC 

30 

Nombre de rencontres entre les 
et C dans la zone profonde. 

animaux 

BN 

ZPBNE 

30 

Nombre de rencontres entre les 
et E dans la zone profonde. 

animaux 

BN 

ZPBNF 

30 

Nombre de rencontres entre les 
et F dans la zone profonde. 

animaux 

BN 

ZPBNG 

30 

Nombre de rencontres entre les 
et G dans la zone profonde. 

animaux 

BN 
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Symbole 

Longueur 

vecteur 

Signification 

ZEBNC 

30 

Nombre de rencontres entre les animaux BN 
et C dans la zone extérieure. 

ZEBNE 

30 

Nombre de rencontres entre les animaux BN 
et E dans la zone extérieure. 

ZEBNF 

30 

Nombre de rencontres entre les animaux BN 
et F dans la zone extérieure. 

ZEBNG 

30 

Nombre de rencontres entre les animaux BN 
et G dans la zone extérieure. 

ZEBSC 

30 

Nombre de rencontres entre les animaux BS 
et C dans la zone extérieure. 

ZEBSE 

30 

Nombre de rencontres entre les animaux BS et 
E dans la zone extérieure. 

ZEBSF 

30 

Nombre de rencontres entre les animaux BS et 
F dans la zone extérieure. 

ZEBSG 

30 

Nombre de rencontres entre les animaux BS et 
G dans la zone extérieue. 

ZPBSC 

30 

Nombre de rencontres entre les animaux BS et 
C dans la zone profonde. 

ZPBSE 

30 

Nombre de rencontres entre les animaux BS et 
E dans la zone profonde 

ZPBSF 

30 

Nombre de rencontres entre les animaux BS et 
F dans la zone profonde. 

ZPBSG 

30 

Nombre de rencontres entre les animaux BS et 
G dans la zone profonde. 

ZEJAN 

30 

Nombre de jeunes du groupe AN dans la zone 
extérieure. 

ZEJAS 

30 

Nombre de jeunes du groupe AS dans la zone 
extérieure. 

ZEJBN 

30 

Nombre de jeunes du groupe BN dans la zone 
extérieure. 

ZEJBS 

30 

Nombre de jeunes du groupe BS dans la zone 


extérieure. 
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Symbole 


Longueur 

vecteur 


Signification 


ZEJJC 

30 

Nombre de jeunes du groupe C dans la 
extérieure. 

zone 

ZEJJE 

30 

Nombre de jeunes du groupe E dans la 
extérieure. 

zone 

ZEJJF 

30 

Nombre de jeunes du groupe F dans la 
extérieure. 

zone 

ZEJJG 

30 

Nombre de jeunes du groupe G dans la 
extérieure. 

zone 

ZPJAN 

30 

Nombre de jeunes du groupe AN dans la 
profonde. 

zone 

ZPJAS 

30 

Nombre de jeunes du groupe AS dans la 
profonde. 

zone 

ZPJBN 

30 

Nombre de jeunes du groupe BN dans la 
profonde. 

zone 

ZPJBS 

30 

Nombre de jeunes du groupe BS dans la 
profonde. 

zone 

ZPJJC 

30 

Nombre de jeunes du groupe C dans la 
profonde. 

zone 

ZPJJE 

30 

Nombre de jeunes du groupe E dans la 
profonde. 

zone 

ZPJJF 

30 

Nombre de jeunes du groupe F dans la 
profonde. 

zone 

ZPJJG 

30 

Nombre de jeunes du groupe G dans la 
profonde. 

zone 

ZET 

30 

Température dans la zone extérieure. 


Z PT 

30 

Température dans la zone profonde. 


ZERH 

30 

Humidité relative dans la zone extérieure. 


ZPRH 

30 

Humidité relative dans la zone profonde. 


ANC 

30 

Nombre de rencontres entre les animaux 
et C. 

AN 

ANE 

30 

Nombre de rencontres entre les animaux 
et E. 

AN 
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Symbole 

Longueur 

vecteur 


Signification 



ANF 

30 

Nombre 
et F. 

de 

rencontres entre 

les 

animaux 

AN 

ANG 

30 

Nombre 
et G. 

de 

rencontres entre 

les 

animaux 

AN 

BNC 

30 

Nombre 
et C. 

de 

rencontres entre 

les 

animaux 

BN 

BNE 

30 

Nombre 
et E. 

de 

rencontres entre 

les 

animaux 

BN 

BNF 

30 

Nombre 
et F. 

de 

rencontres entre 

les 

animaux 

BN 

BNG 

30 

Nombre 
et G. 

de 

rencontres entre 

les 

animaux 

BN 

BSC 

30 

Nombre 
et C. 

de 

rencontres entre 

les 

animaux 

BS 

BSE 

30 

Nombre 
et E. 

de 

rencontres entre 

les 

animaux 

BS 

BSF 

30 

Nombre 
et F. 

de 

rencontres entre 

les 

animaux 

BS 

BSG 

30 

Nombre 
et G. 

de 

rencontres entre 

les 

animaux 

BS 

PRMGAN 

30 

Nombre de proies mangées 
AN, en une unité de temps. 

par 

les animaux 

PRMGBN 

30 

Nombre de proies mangées 
BN, en une unité de temps. 

par 

les animaux 

PRMGBS 

30 

Nombre 
BS, en 

de proies mangées 
une unité de temps. 

par 

les animaux 

PRMGAS 

30 

Nombre 
AS, en 

de proies mangées 
i une unité de temps. 

par 

les animaux 

AN 

30 

Nombre 

total d'animaux AN. 




AS 

30 

Nombre 

total d'animaux AS. 




BN 

30 

Nombre 

total d'animaux BN. 




BS 

30 

Nombre 

total d'animaux BS. 






64 


R. R. TERCAFS. — SIMULATION DUN SYSTÈME ÉCOLOGIQUE 


47, 19 


Symbole 

Longueur 

vecteur 

Signification 

C 

30 

Nombre total d’animaux C. 

E 

30 

Nombre total d’animaux E. 

F 

30 

Nombre total d'animaux F 

G 

30 

Nombre total d'animaux G. 

T 

30 

Indice temps, utilisé dans les graphiques unique¬ 
ment. 

TAS 

30 

Variable intermédiaire intervenant dans 
VARIAT. 

TBS 

30 

Variable intermédiaire intervenant dans 
VARIAT. 

GGAMJA 

30 

Nombre de jeunes du groupe A nécessaire pour 
obtenir une proie normale. 

GGAMJB 

30 

Nombre de jeunes du groupe B nécessaire pour 
constituer une proie normale. 

GGAMJC 

30 

Nombre de jeunes du groupe C nécessaire pour 
constituer une proie normale. 

GGAMJE 

30 

Nombre de jeunes du groupe E nécessaire pour 
constituer une proie normale. 

GGAMJF 

30 

Nombre de jeunes du groupe F nécessaire pour 
constituer une proie normale. 

GGAMJG 

30 

Nombre de jeunes du groupe G nécessaire pour 
constituer une proie normale. 

VVJA 

30 

Vitesse de déplacement des jeunes du groupe A. 

VVJB 

30 

Vitesse de déplcaement des jeunes du groupe B. 

VVJC 

30 

Vitesse de déplacement des jeunes du groupe C. 

VVJE 

30 

Vitesse de déplacement des jeunes du groupe E. 

WJ F 

30 

Vitesse de déplacement des jeunes du groupe F. 

WJ G 

30 

Vitesse de déplacement des jeunes du groupe G. 

EEPSJA 

30 

Coefficient de correction prédateur-proie pour 


les jeunes du groupe A. 
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Symbole 

Longueur 

vecteur 

Signification 

EEPSJB 

30 

Coefficient de correction prédateur-proie pour 
les animaux du groupe B. 

EEPSJC 

30 

Variable non utilisée dans le présent programme. 

DDJA 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe A. 

DDJB 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe B. 

DDJC 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe C. 

DDJE 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe E. 

DDJF 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe F. 

DDJG 

30 

Diamètre de perception ou d’influence des jeunes 
du groupe G. 

J AN 

30 

Nombre total de jeunes du groupe AN. 

JAS 

30 

Nombre total de jeunes du groupe AS. 

J BS 

30 

Nombre total de jeunes du groupe BS. 

JBN 

30 

Nombre total de jeunes du groupe BN. 

lie 

30 

Nombre total de jeunes du groupe C. 

UE 

30 

Nombre total de jeunes du groupe E. 

JJF 

30 

Nombre total de jeunes du groupe F. 

JJG 

30 

Nombre total de jeunes du groupe G. 

EGMBNH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe BN. 

EGMBNT 

100 

Gamma thermique'dans la zone extérieure pour 
le groupe BN. 

ABBNT 

100 

Abscisses du thermopréférendum pour le groupe 

BN. : - .3-. 

RPBNT 

100 

■Coefficient de répartition-température pour le 
groupe BN. 
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Symbole 

Longueur 

vecteur 

Signification 

RPBNH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe BN. 

ABSBNT 

100 

Abscisses du gamma thermique pour le groupe 
BN. 

ABSBNH 

100 

Abscisses du gamma hygrométrique pour le 
groupe BN. 

EGBNTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe BN. 

EGBNHG 

100 

Gamma hygrométrique globalisé dans la zone 
extérieure pour le groupe BN. 

PGBNTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe BN. 

PGBNHG 

100 

Gamma hygrométrique globalisé dans la zone 
profonde pour le groupe BN. 

EGMANH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe AN. 

EGMANT 

100 

Gamma thermique dans la zone extérieure pour 
le groupe AN. 

ABANT 

100 

Abscisses du thermopréférendum pour le groupe 
AN. 

RPANT 

100 

Coefficient de répartition-température pour le 
groupe AN. 

RPANH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe AN. 

ABSANT 

100 

Abscisses du gamma thermique pour le groupe 
AN. 

ABSANH 

100 

Abscisses du gamma hygrométrique pour le 
groupe AN. 

EGANTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe AN. 

EGANHG 

100 

Gamma hygrométrique globalisé dans la zone 


extérieure pour le groupe AN. 
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Symbole 

Longueur 

vecteur 

Signification 

PGANTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde du groupe AN. 

PGANHG 

100 

Gamma hygrométrique globalisé pour le groupe 
AN dans la zone profonde. 

EGMASH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe AS. 

EGMAST 

100 

Gamma thermique dans la zone extérieure pour 
le groupe AS. 

ABAST 

100 

Abscisses du thermopréférendum pour le groupe 

AS. 

RPAST 

100 

Coefficient de répartition-température pour le 
groupe AS. 

RPASH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe AS. 

ABSAST 

100 

Abscisses du gamma thermique pour le groupe 

AS. 

RPAST 

100 

Coefficient de répartition-température pour le 
groupe AS. 

RPASH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe AS. 

ABSAST 

100 

Abscisses du gamma thermique pour le 
groupe AS. 

ABSASH 

100 

Abscisses du gamma hygrométrique pour le 
groupe AS. 

EGASTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe AS. 

EGASHG 

100 

Gamma hygrométrique globalisé dans la zone 
extérieure pour le groupe AS. 

PGASTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe AS. 

PGASHG 

100 

Gamma hygrométrique globalisé dans la zone 


profonde pour le groupe AS. 
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Symbole 

Longueui 

vecteur 

Signification 

EGMBSH 

100 

Gamma hygrométrique pour la zone extérieure 
pour le groupe BS. 

EGMBST 

100 

Gamma thermique globalisé pour la zone exté¬ 
rieure pour le groupe BS. 

ABBST 

100 

Abscisses du thermopréférendum pour le groupe 

BS. 

RPBST 

100 

Coefficient de répartition-température pour le 
groupe BS. 

RPBSH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe BS. 

ABSBST 

100 

Abscisses du gamma thermique pour le groupe 

BS. 

ABSBSH 

100 

Abscisses du gamma hygrométrique pour le 
groupe BS. 

EGBSTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe BS. 

EGBSHG 

100 

Gamma hygrométrique globalisé dans la zone 
extérieure pour le groupe BS. 

PGBSTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe BS. 

PGBSHG 

100 

Gamma hygrométrique globalisé pour la zone 
profonde pour le groupe BS. 

EGMCCH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe C. 

EGMEEH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe E. 

EGMFFH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe F. 

EGMGGH 

100 

Gamma hygrométrique dans la zone extérieure 
pour le groupe G. 

EGMCCT 

100 

Gamma thermique dans la zone extérieure pour 
le groupe C. 

EGMEET 

100 

Gamma thermique dans la zone extérieure pour 


la zone C. 
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Symbole 

Longueur 

vecteur 

Signification 

EGMFFT 

100 

Gamma thermique dans la zone extérieure pour 
la zone E. 

EGMGGT 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe G. 

ABCCT 

100 

Abscisses du thermopréférendum du groupe C. 

ABEET 

100 

Abscisses du thermopréférendum du groupe E. 

ABFFT 

100 

Abscisses du thermopréférendum du groupe F. 

ABGGT 

100 

Abscisses du thermopréférendum du groupe G. 

RPCCT 

100 

Coefficient de répartition-température pour le 
pour le groupe C. 

RPEET 

100 

Coefficient de répartition-température pour le 
groupe E. 

RPFFT 

100 

Coefficient de répartition-température pour le 
groupe F. 

RPGGT 

100 

Coefficient de répartition-température pour le 
pour le groupe G. 

RPCCH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe C. 

RPEEH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe E. 

RPFFH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe F. 

RPGGH 

100 

Coefficient de répartition-humidité pour le grou¬ 
pe G. 

ABSCCT 

100 

Abscisses du gamma thermique pour le groupe 
C. 

ABSEET 

100 

Abscisses du gamma thermique pour le qroupe 
E, 

ABSFFT 

100 

Abscisses du gamma thermique pour le groupe 
F. 

ABSGGT 

100 

Abscisses du gamma thermique pour le groupe 
G. 
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Symbole 

Longueur 

vecteur 

Signification 

ABSCCH 

100 

Abscisses du gamma hygrométrique pour le 
groupe C. 

ABSEEH 

100 

Abscisses du gamma hygrométrique pour le 
groupe E. 

ABSFFH 

100 

Abscisses du gamma hygrométrique pour le 
groupe F. 

ABSGGH 

100 

Abscisses du gamma hygrométrique pour le 
groupe G. 

EGCCTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe C. 

EGEETG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe E. 

EGFFTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe F. 

EGGGTG 

100 

Gamma thermique globalisé dans la zone exté¬ 
rieure pour le groupe G. 

EGCCHG 

100 

Gamma hygrométrique globalisée dans la zone 
extérieure pour le groupe C. 

EGEEHG 

100 

Gamma hygrométrique globalisée dans la zone 
extérieure pour le groupe E. 

EGFFHG 

100 

Gamma hygrométrique globalisée dans la zone 
extérieure pour le groupe F. 

EGGGHG 

100 

Gamma hygrométrique globalisée dans la zone 
extérieure pour le groupe G. 

PGCCTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe C. 

PGEETG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe E. 

PGFFTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe F. 

PGGGTG 

100 

Gamma thermique globalisé dans la zone pro¬ 
fonde pour le groupe G. 

PGCCHG 

100 

Gamma hygrométrique globalisé dans la zone 


profonde du groupe C. 
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Symbole 

Longueur 

vecteur 

Signification 

PGEEHG 

100 

Gamma hygrométrique globalisé dans la zone 
profonde du groupe E. 

PGFFHG 

100 

Gamma hygrométrique globalisé dans la zone 
profonde du groupe F. 

PGGGHG 

100 

Gamma hygrométrique globalisé dans la zone 
profonde du groupe G. 

RTET 

100 

Coefficient de répartition dû à la température 
(zone ext.). 

RTPT 

100 

Coefficient de répartition dû à la température 
(zone profonde). 

RTEH 

100 

Coefficient de répartition dû à l'humidité rela¬ 
tive (zone ext.). 

RTPH 

100 

Coefficient de répartition dû à l'humidité rela¬ 
tive (zone profonde). 

MAN 

1 

Coefficient de mortalité affectant le groupe AN. 

MAS 

1 

Coefficient de mortalité affectant le groupe AS. 

MBN 

1 

Coefficient de mortalité affectant le groupe BN. 

MBS 

1 

Coefficient de mortalité affectant le groupe BS. 

MC 

1 

Coefficient de mortalité affectant le groupe C. 

ME 

1 

Coefficient de mortalité affectant le groupe E. 

MF 

1 

Coefficient de mortalité affectant le groupe F. 

MG 

1 

Coefficient de mortalité affectant le groupe G. 


III. RESULTATS 

L'analyse détaillée des résultats obtenus par cette méthode de simu¬ 
lation sera donnée dans la deuxième partie de ce travail. Néanmoins, nous 
présentons dans les figures 1 et 2 un exemple de tracé. 

Nous avons simulé un milieu écologique constitué de 2 zones (zone 
extérieure et profonde) où les conditions climatiques sont différentes et 
varient en fonction du temps. Suivant leurs préférendums, les animaux 
choisissent la zone leur convenant le mieux. 

La période expérimentale est divisée en 30 unités de temps. Les 
valeurs des données sont reprises dans le table 1. 



72 


R. R. TERCAFS. — SIMULATION D UN SYSTÈME ÉCOLOGIQUE 


47, 19 


REMERCIEMENTS 


Nous remercions M. le Professeur Paul Berthet (Université de 
Louvain) d’avoir bien voulu nous conseiller en divers points de ce 
travail. 

Les passages en ordinateur ont été assurés par le Centre de Calcul 
de l'Université de Liège. 


TABLEAU 1 

Données utilisées pour la simulation d'un système écologique à 2 zones climatiquement 
différentes (pour la signification des abréviations, voir pages 47 et suivantes). 
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RESUME 

Un modèle écologique simulant les relations entre 2 populations de 
carnivores et 4 populations d'herbivores vivant dans 2 zones climatique¬ 
ment différentes a été développé. 

Après discussion des modèles mathématiques précédemment décrits, 
on détaille les paramètres physiques et biologiques utilisés, ainsi que leurs 
relations. Ces relations sont quantifiées sous forme de plusieurs systèmes 
d’équations. Les systèmes sont programmés en FORTRAN IV. 35 sous- 
programmes sont utilisés : leur rôle et le détail de leur structure est donné, 
ainsi que les procédures d’assemblage pour les ordinateurs IBM 360/65 
et 360/44. 

Des exemples des résultats obtenus sont donnés. 

SUMMARY 

An ecological model simulating the relations between 2 carnivorous 
and 4 herbivorous populations living in 2 areas of different climate was 
developed. 
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After discussion of mathematical models previously described, one 
details the physical and biological parameters used, and also their 
relations. These relations are quantified in some équation Systems. 
Systems are programmed in FORTRAN IV. 

35 subroutines are used : their object and their structure is given, with 
the assembly procedure for IBM 360/65 and 360/44 computers. 

Some examples of results obtained by these methods are given. 

Institut royal des Sciences naturelles de Belgique. 
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